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. Context. I Zw 18, ever since regarded as the prototypical blue compact dwarf (BCD) galaxy, is, quite ironically, the most atypical 

BCD known. This is because its large low-surface brightness (LSB) envelope is not due to an old underlying stellar host, as invariably 
is the case for typical BCDs, but entirely due to extended nebular emission (Papaderos et al. 2002; hereafter P02). 
Aims. Our goal is to explore I Zw 18 and its detached C component I Zw 18 C down to an unprecedently faint surface brightness p 
(mag/d") level in order to gain further insight into the structural properties and evolutionary history of this enigmatic galaxy pair. 
Methods. We present a photometric analysis of the entire set of archival HST ACS V, R and / band data for I Zw 18. 
Results. Radial color profiles for I Zw 18 C reveal blue and practically constant colors (0±0.05) down to p~21 .6, and a previously 
[ undisclosed, slightly redder (V-I~0.2), stellar population in its extreme periphery (p ~29). We argue that stellar diffusion over t~10 8 

O j. yr and the associated stellar mass filtering effect (P02) can consistently account for the observed properties of the stellar component 

I ■ in the outskirts of I Zw 18 C. This process, in combination with propagating star formation with a mean velocity of ~20 km s" 1 can 

reproduce all essential characteristics of I Zw 18 C within ~ t. An extremely faint substrate of older stars can neither be ruled out nor 
does need be postulated. As for I Zw 18, we find that nebular emission (ne) extends out to ~16 stellar scale lengths, shows a nearly 
exponential outer profile, and provides at least 1/3 of the total optical emission, ne dominates already at p ~23.5, as evident from e.g. 
the uniform and extremely blue (V-/«-l, /J-/«-1.4) colors of the LSB envelope of I Zw 18. 

Conclusions. The case of I Zw 18 suggests caution in studies of distant galaxies in dominant stages of their evolution, rapidly 
assembling their stellar mass at high specific star formation rates (SSFRs). It calls attention to the fact that ne is not necessarily 
cospatial with the underlying ionizing and non-ionizing stellar background, neither has to scale with its surface density. The prodigious 
energetic output during dominant phases of galaxy evolution may result in large exponential ne envelopes, extending much beyond 
the still compact stellar component, just like in I Zw 18. Therefore, the morphological paradigm of I Zw 18, while probably unique 
in the nearby Universe, may be ubiquitous among high-SSFR galaxies at high redshift. Using I Zw 18 as reference, we show that 
extended ne may introduce substantial observational biases and affect several of the commonly studied fundamental galaxy relations. 
Among others, we show that the surface brightness profiles of distant morphological analogs to I Zw 18 may be barely distinguishable 
' from Sersic profiles with an exponent 2S r\ 5=5, thus mimicking the profiles of massive galaxy spheroids. 

Key words, galaxies: dwarf - galaxies: starburst - galaxies: structure - galaxies: evolution - galaxies: evolution - galaxies: high- 
redshift 
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1 . Introduction the unconquered prototypical example of this enigmatic galaxy 

class. 



C3 Even four decades after its discovery (Sargent & Searle, 1970), 

the blue compact dwarf (BCD) galaxy I Zw 18 continues to I Zw 18 was originally described by Zwicky (1966) as a 

attract considerable interest and feed intense debates in extra- pair of compact galaxies, later on recognized to be SF regions 

galactic research. Its low oxygen abundance (Searle & Sargent, within the same galaxy, the brighter northwestern (NW) and the 

1972), established in numerous subsequent studies (Lequeux et fainter southeastern (SE) component, separated by ^6" (cf Fig. 

al., 1979; French, 1980; Kinman & Davidson, 1981; Pagel et 1). Subsequent work has shown that these regions are embed- 

al., 1992; Skillman & Kennicutt, 1993; Martin, 1996; Vflchez ded within an extended, low-surface brightness (LSB) envelope 

& Iglesias-Paramo, 1998; Izotov & Thuan , 1998a,b; Izotov et (Davidson et al., 1989; Dufour & Hester, 1990; Ostlin et al., 

al., 1999, among others) to be 12+log(0/H)«7.2, makes it the 1996; Dufour et al., 1996a; Martin, 1996), whose rich filamen- 

third most metal-poor star-forming (SF) galaxy in the nearby tary substructure was impressively revealed with the advent of 

Universe, after SBS 0335-052 W (Izotov et al., 2005; Papaderos the HST (Hunter & Thronson, 1995; Dufour et al., 1996b). That 

et al., 2006a; Izotov et al., 2009) and DDO68 (Pustilnik et al., nebular emission (hereafter ne) is very strong in the central part 

2005; Izotov & Thuan, 2007). Despite a meanwhile long record of I Zw 18 and its north-western super-shell was spectroscopi- 

of extremely metal-poor (12+log(0/H)^7.6) BCDs (hereafter cally documented early on. For example, Izotov et al. (2001a) 

XBCDs) discovered in the recent years (see e.g. Papaderos et have shown on the basis of deep long slit spectroscopy that the 

al., 2008; Guseva et al., 2009, for a review), I Zw 18 remains equivalent width (EW) of the Ha emission line rises to 1700 A 

northwest of region NW and that ne is present as far away as 15" 

from it (regions labeled "Har arc" and "Loop" in Fig. 1). The 

Send offprint requests to: P. Papaderos; papaderos@astro.up.pt EW(Hd') morphology of I Zw 18 was first studied with high- 
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resolution ground-based imagery by Ostlin et al. (1996) who de- 
scribed a horseshoe-shaped rim of intense (EW(Hff)-1500 A) 
ne encompassing region NW. This conspicuous EW pattern was 
later on confirmed through HST WFPC2 data (Papaderos et al., 
2001; Izotov et al., 2001a) and, more impressively, by Vflchez & 
Iglesias-Paramo (1998) who were the first to present a 2D spec- 
troscopic study of the chemical abundance patterns of the warm 
interstellar medium (ISM) in I Zw 18. 

Much less attention has been attracted by the fainter de- 
tached C component of I Zw 18 (hereafter I Zw 18C), located 
~ 22" northwest of region NW. Dufour et al. (1996a), Petrosian 
et al. (1997), Izotov & Thuan (1998a), van Zee et al. (1998) 
and Izotov et al. (2001a) have shown it to have the same re- 
cession velocity as the main body, thus establishing its physical 
association to I Zw 18. This was also shown through interfero- 
metric 21cm studies (van Zee et al., 1998, see also Viallefond 
et al. (1987)) which revealed that I Zw 18 and I Zw 18 C are 
immersed within a large common HI complex with a projected 
size of 60" x 45" connecting with a southern tail with no 
optical counterpart. The SF activity in I Zw 18 C is known to 
be weak with its EW(Ho') not exceeding ~60 A along its major 
axis (Izotov et al., 2001a, see also van Zee et al. (1998)). Despite 
deep Keck II spectroscopy, Izotov et al. (2001a) failed to detect 
oxygen lines, so its oxygen abundance is not known. 

Traditionally, the distance to I Zw 18 has been taken to 
be 10 Mpc, assuming a pure Hubble flow recessional veloc- 
ity. However, Izotov et al. (1999) have argued based on an 
HST color-magnitude diagram (CMD) study that the distance 
to I Zw 18 has to be at least 15 Mpc, and most likely ~20 Mpc, 
in order for its brightest stars being massive enough to account 
for the ionizing flux observed. This upper distance value has re- 
cently received independent support by Fiorentino et al. (2010). 
These authors have identified three long-period Cepheid candi- 
dates in I Zw 18, which, if interpreted as classical Cepheids, 
imply by the Wesenheit relation a distance of 19.0^2 j Myr. In 
the following, we shall throughout adopt a distance D=19 Mpc 
to both I Zw 1 8 and I Zw 1 8 C and convert distance-dependent 
quantities from the literature accordingly. It should be noted, 
however, that the assumed distance has practically no influence 
on the main conclusions from this study. 

The wealth of dedicated studies of I Zw 18 highlight the im- 
portance placed on this XBCD as precious nearby laboratory for 
exploring collective star formation and the associated feedback 
process under metallicity conditions approaching those in distant 
protogalactic systems. Some examples include the consideration 
of I Zw 18 as reference object for many dwarf galaxy chemical 
evolution models (Matteuchi & Tosi, 1985; Roy & Kunth, 1995; 
Mas-Hesse & Kunth, 1999; Legrand, 2000; Recchi et al., 2004, 
among others), the great deal of effort put in the determination of 
its chemical abundance patterns in its neutral ISM (Kunth et al., 
1994; Aloisi et al., 2003; Lecavelier des Etangs et al., 2004), and 
in the study of its dust and molecular gas content (e.g. Cannon 
et al., 2002; Leroy et al., 2007), the thorough exploration of the 
excitation mechanisms of its brightest Hll region (Stasihska & 
Schaerer, 1999; Pequignot, 2008), and the deep spectroscopic 
studies that led to the discovery of Wolf-Rayet stellar features in 
it (Izotov et al., 1997b; Legrand et al., 1997). 

But arguably, the most longstanding debate associated with 
I Zw 18 ever since its discovery concerns its evolutionary sta- 
tus. In this regard, various interpretations have been put for- 
ward, ranging from I Zw 18 being a bona fide young galaxy, 
currently forming its first stellar generation (Sargent & Searle, 
1970, see also Izotov & Thuan (1999)), all through the diamet- 



rically opposite picture of an ancient "slowly cooking" dwarf 
galaxy that is forming stars continuously over the Hubble time 
(Legrand, 2000; Legrand et al., 2001). Notwithstanding an im- 
pressive amount of high-quality multiwavelength data and nu- 
merous dedicated analyses of considerable effort and sophistica- 
tion, the convergence towards a consensual view on the evolu- 
tionary status of I Zw 18 has been slow. 

CMD analyses, based on HST data, have primarily been fo- 
cusing on the question of whether or not I Zw 18 contains a size- 
able population of evolved red giant branch (RGB) stars, simi- 
lar to typical (12+log(0/H)>8, see e.g. Kunth & Ostlin, 2000) 
BCDs. In the latter, an extended envelope of resolved RGB stars 
around the SF component (e.g. Tosi et al., 2001; Crone et al., 
2002) nicely echoes the since-long observationally established 
fact of an evolved underlying host galaxy in these systems (e.g. 
Loose & Thuan, 1986; Papaderos et al., 1996a; Cairos et al., 
2001; Bergvall & Ostlin, 2002; Noeske et al., 2003; Gil de Paz 
& Madore, 2005). Initial CMD studies suggested for the main 
body an age between several 10 Myr (Hunter & Thronson, 1995; 
Dufour et al., 1996b) and ~ 1 Gyr (Aloisi et al., 1999). Ostlin 
(2000) argued from an HST NICMOS near infrared (NIR) study 
that a fit to the J vs. / - H CMD of I Zw 18 is best achieved 
for a stellar population of age ~5 Gyr. The subsequent identifi- 
cation of five carbon star candidates with an estimated age 0.5-1 
Gyr by Ostlin & Mouhcine (2005) is in accord with that con- 
clusion, even though the number of evolved star candidates in 
all above studies (about a dozen altogether) was recognized to 
be surprisingly small compared to typical BCDs of equal lumi- 
nosity. A significant step forward has been possible through the 
advent of HST ACS, allowing (Izotov & Thuan, 2004, hereafter 
IT04) to extend point source photometry to magnitudes as faint 
as 29 mag in the V and / band and revisit the question of the 
presence of RGB stars in I Zw 18. IT04 found, in addition to nu- 
merous blue main sequence and blue and red supergiants with an 
age S5100 Myr, an older population of asymptotic giant branch 
(AGB) stars with an age between 0.1 and 0.5 Gyr. This study, 
in which no RGB were detected, has been the first to also ex- 
plore the spatial distribution of stars of different ages in I Zw 18. 
The upper age limit of 0.5 Gyr for the oldest stars in I Zw 18 
(IT04) was subsequently relaxed from a re-analysis of the same 
data by Yakobchuk & Izotov (2006) which revealed an untyp- 
ically small number of RGB candidates. Various other efforts 
have been made to improve on the CMD analysis of IT04 by 
pushing point source photometry to by 1-2 mag fainter levels 
(Momany et al., 2005; Tosi et al., 2007; Tikhonov, 2007). The 
faintest (>29 mag) point sources in those CMDs cover almost 
uniformly the color range between <-l mag and >2 mag. It is 
worth pointing out that, whereas divergent in their conclusions 
regarding stellar age, all CMD analyses for I Zw 18 consistently 
indicate a conspicuous absence of an extended stellar LSB en- 
velope surrounding regions NW&SE, at sharp contrast to any 
previously studied BCD. 

As for I Zw 18 C, CMD analyses yield an upper age between 
a few ten and hundred Myr (Dufour et al., 1996b; Aloisi et al., 
1999). Recently, Jamet et al. (2010) employing a probabilistic 
CMD modeling technique reported an upper age of ~125 Myr, 
without, however, strictly ruling out the presence of older stars. 
An age of the same order was previously inferred for I Zw 18 C 
from a combined CMD and evolutionary spectral synthesis study 
by Izotov et al. (2001a). 

From the viewpoint of surface photometry, diametrically dif- 
ferent conclusions on the photometric structure and evolutionary 
status of I Zw 1 8 were drawn by Kunth & Ostlin (2000, hereafter 
KOOO) and Papaderos et al. (2002, hereafter P02). Nevertheless, 
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Fig. 1. Three-color composite image of I Zw 18 and I Zw 18 C, combining HST ACS data in V, R and / (blue, green and red channel, respectively). 
The position of the star-forming regions NW and SE is indicated by crosses. The regions labeled "Loop" and "Ha arc" were studied through 
deep Keck II long slit spectroscopy by Izotov et al. (2001a). The region cj at the southeastern tip of I Zw 18 shows comparatively weak nebular 
emission, its colors allow therefore to place meaningful constraints on the age of the stellar component (Papaderos et al., 2002). In the magnified 
version of region NW (inset), combining the unsharp masked images I c , R c and V„ (see Sect. 3 for details), about 30 point sources, surrounded by 
a complex network of ionized gas shells are discernible. The irregular blue galaxy I Zw 18 C is located ~22" northwest of region NW (x2 kpc at 
the assumed distance of 19 Mpc to I Zw 18). It shows faint nebular emission in its bluer southeastern tip and central star cluster complex C (Izotov 
et al., 2001a). North is at the top and east to the left. 



these two studies were the first to demonstrate on the basis of 
surface photometry that I Zw 18 is not presently forming its first 
stellar generation but contains a substantial unresolved stellar 
background of intermediate age. 

As a matter of fact, much of the disparity between these stud- 
ies has been due to the different importance they ascribed to the 
presence and photometric impact of ne. KO00 concluded that SF 
activity in 1 Zw 18 is hosted by an old, extended stellar disk that 
dominates the stellar mass, just like in typical BCDs. Their ratio- 
nale has mainly been based on their finding that I Zw 18 shows 
an exponential intensity decrease and reddish (B-R^Q.6 mag) 
colors in its LSB envelope (9"^i?*^20")- On the assumption 



that stellar emission dominates throughout, this color translates 
by a continuous star formation model to an age of ^5 Gyr. The 
central surface brightness /iq and exponential scale length a of 
the disk, read off the B band surface brightness profile (SBP) of 
KO00, imply that the old stellar host contains ~l/2 of the emis- 
sion and the bulk of the stellar mass in I Zw 18. Note that the 
stellar disk interpretation for I Zw 18 is in qualitative agreement 
with the evolutionary scenario by Legrand et al. (2001, see also 
Legrand (2000)). These authors argued that the low and uniform 
gas-phase metallicity of I Zw 18 is reproducible through contin- 
uous low-level star formation throughout the main HI complex 
of I Zw 18 (45"x60") over the past -14 Gyr. This process would 
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produce an extended stellar disk of extremely low surface bright- 
ness (Jl ^28 B mag/D"). 

P02 (see also Papaderos et al., 2001), on the other hand, 
called into question the conclusions by KO00 by invoking var- 
ious lines of observational evidence. First, they have empiri- 
cally shown that an exponential outer intensity drop off is a 
generic property of the nebular halo of starbursting dwarf galax- 
ies. Consequently, the exponentiality of the LSB envelope of 
I Zw 18 is not per se a compelling argument for it being due 
to a stellar disk. This also applies to its reddish colors which 
can naturally be accounted for by photoionized gas of subsolar 
metallicity (see e.g. Kriiger et al., 1995). Secondly, and in a more 
straight forward approach, P02 used HST WFPC2 narrow band 
images to bidimensionally subtract the [Olll] and Ha line emis- 
sion from broad band HST data in order to isolate and study 
the residual stellar emission in I Zw 18. This correction led to 
the virtual removal of the filamentary LSB envelope, proving its 
gaseous nature. Specifically, P02 have shown that ne dominates 
the line-of-sight intensity already at a photometric radius R* ^6" 
and contributes between 30% and 50% of the R band luminosity 
of I Zw 18. Broad band images, after decontamination from neb- 
ular line emission (though still affected by nebular continuum 
emission), were then used to study the photometric structure and 
color distribution of the stellar component of I Zw 18. SBPs 
computed from them have revealed a very compact host which, 
at sharp contrast to typical BCDs, shows practically no radial 
color gradients and overall very blue (^0.1 mag) colors down 
to a limiting surface brightness fi ~26 mag/D". These excep- 
tional properties were interpreted as evidence for youth: young 
stars in I Zw 18 did not have had enough time to migrate signifi- 
cantly far from heir initial locus and gradually form the extended 
stellar host that is typical of evolved BCDs. This, and the blue 
optical and NIR colors of the southeastern tip of I Zw 1 8 (region 
to in the notation of P02, cf Fig. 1) where ne is weak led P02 
to conclude that most of the stellar mass in I Zw 18 has formed 
within the past 0.5 Gyr. Hence, the picture put forward by P02 
is that I Zw 18 is a cosmologically young object that presently 
undergoes its dominant formation phase and contains a small, if 
any, mass fraction of stars older than ~1 Gyr. Further support to 
this conclusion came from a subsequent NIR study of I Zw 1 8 by 
Hunt et al. (2003). As for I Zw 18 C, the nearly constant blue col- 
ors («0 mag) determined within its Holmberg radius, lent further 
support to the youth interpretation that was previously advocated 
by Izotov et al. (2001a). 

However, important aspects of the I Zw 18 system could not 
be conclusively addressed from previous photometric studies. 
For example, whereas SBPs for I Zw 18 C by P02 reach a surface 
brightness level p. ~ 28 B mag/D", their large photometric uncer- 
tainties already below *26 mag/D" have practically prevented an 
assessment of the question of whether a redder underlying stel- 
lar population dominates in the extremely faint periphery of the 
galaxy. Clearly, this issue is central to the understanding of the 
evolutionary status of I Zw 18 C. One may argue that, since the 
evolved stellar host generally dominates for p ^ 24.5 B mag/D" 
(Papaderos et al., 1996a, hereafter P96a), it would have been de- 
tected in I Zw 18 C, should have been present. This is a circular 
argument, however, given that empirical relations established for 
typical BCDs should not be taken for granted for young dwarf 
galaxy candidates. Similarly, due to the shallowness of previous 
surface photometry no definite conclusions could be drawn re- 
garding the ultra-LSB disk predicted by Legrand (2000). 

As for the main body of I Zw 18, previous studies did not 
had the sensitivity to pin down the maximal extent, morphol- 
ogy and color pattern of the LSB envelope, adding potentially 



important constraints to chemodynamical and spectrophotomet- 
ric models for I Zw 18. From such considerations, extremely 
deep surface photometry appears crucially important for further 
advancing our understanding of the photometric structure and 
evolutionary status of I Zw 18 and I Zw 18 C. This is particu- 
larly true for deep / band surface photometry which is entirely 
lacking both for I Zw 18 and I Zw 18 C. Note that / band HST 
WFPC2 images included the main body of I Zw 18 only and 
were not deep enough for a study of its LSB envelope. This data 
set was therefore not used in the surface photometry analysis by 
P02. 

This study is motivated by the availability of an unprece- 
dently deep set of archival HST ACS V, R and / broad band data 
that has accumulated over the past few years. / band photometry 
is an important asset in this respect, not only due to its higher 
sensitivity to a putative old stellar background but also because 
it offers, in combination with V and R data, a sensitive tracer of 
ne both in the center and the LSB envelope of I Zw 18. This is 
because the / band is affected by nebular continuum emission 
only, whereas the V and R transmission curves additionally in- 
clude, respectively, the strong [Olll]/l/l4959,5007 and Ha emis- 
sion lines. As a result, regions with strongly enhanced ne can 
readily be identified by their extremely blue (-0.5 . . .-1.4) V-I 
and R-I colors (see e.g. Papaderos et al., 1998, P02 and refer- 
ences therein for a discussion and examples among XBCDs). 

This paper is organized as follows: in Sect. 2 we discuss 
the data processing and SBP derivation technique used and in 
Sect. 3 the structural and morphological properties of I Zw 18 C 
and I Zw 18. Section 4 concentrates on the evolutionary status 
and the formation process of I Zw 18 C (Sect. 4.1) under con- 
sideration of the effect that the diffusion of young stars would 
have on the observed colors (Sect. 4.1.1). The evolutionary sta- 
tus of I Zw 18 and the hypothesis of an ultra-LSB underlying 
stellar disk are discussed on the basis of the present photometric 
analysis in Sects. 4.2 and 4.2.1, respectively. In Sect. 4.3 we use 
I Zw 18 as template to briefly explore the biases that extended ne 
may introduce in photometric studies of morphologically anal- 
ogous star-forming galaxies at higher redshift (z). The main re- 
sults from this study are summarized in Sect. 5. 

2. Data processing 

This study is based on the entire set of archival HST ACS broad 
band images for I Zw 18 that has been acquired through the HST 
programs 9400 (PI: Thuan) and 10586 (PI: Aloisi). It comprises 
38, 65 and 81 images in the filters F555W (V), F606W (broad 
VR, referred to in the following as R) and F814W (I), sum- 
ming up to on-source exposures of 87, 55 and 101 ksec, respec- 
tively. This is the deepest imaging data set currently available for 
I Zw 18, with an integration time in R and / equaling ~ 1/3 of the 
time spent on the HST ACS Ultra Deep Field (Beckwith et al., 
2006) in the filters F606W and F775W. 

The data processing was carried out using IRAF 1 and ESO 
MIDAS 2 . Photometric quantities refer to the Vega system. 

Since the main goal of this study is deep surface photome- 
try, its most critical aspect is the removal of diffuse and com- 
pact background sources (diffuse extragalactic background, zo- 

1 IRAF is the Image Reduction and Analysis Facility distributed by 
the National Optical Astronomy Observatory, which is operated by the 
Association of Universities for Research in Astronomy (AURA) under 
cooperative agreement with the National Science Foundation (NSF). 

2 Munich Image Data Analysis System, provided by the European 
Southern Observatory (ESO) 
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Fig. 2. Combined / band exposure of I Zw 18 and I Zw 18 C with point- 
and compact sources in the vicinity of the galaxies (polygonal regions 
A and B) marked with rectangles. North is up and east to the left. 



diacal light, foreground stars and background galaxies, respec- 
tively). After subtraction of the diffuse background, we therefore 
checked to which extent compact or point sources (cps) within 
the extended LSB envelope of 1 Zw 1 8 can affect SBPs and color 
profiles. For example, for an extended source of constant sur- 
face brightness fi = 29 mag/D", the total apparent magnitude 
m of a circular annulus with 19"<i?*<20" (roughly the radius 
of I Zw 18) is 23.8 mag. At this intensity level, already a sin- 
gle faint (m=25 mag) background cps can introduce an error of 
0.3 mag in surface photometry. We adopted the following proce- 
dure: after image alignment and correction for cosmics, we used 
the combined images in the three filters to compile a catalogue of 
cps in the relevant portion of the field of view. In doing this, we 
disregarded cps in I Zw 18 (roughly the area subtended by the 25 
R' mag/D" isophote in Fig. 8) and in I Zw 18 C, as well as com- 
pact clumps of ne, identified by their blue V-I and R-I color. All 
cps in each frame were in turn replaced by the mean intensity in 
the adjacent area. Special care was given to the removal of two 
background galaxies close to the western super-shell of I Zw 18. 
This was done by subtracting a 2D model, computed with the 
Bender & Mollenhoff (1987) algorithm. In total, -140 cps were 
removed in the field of interest around I Zw 1 8 and I Zw 1 8 C 
(polygonal regions labeled A and B in Fig. 2). Their integral 7 
magnitude of 19.17 mag and 21.3 mag within the regions consid- 
ered (~1600 □" and ~200 □" for A and B, respectively) corre- 
sponds to a mean surface brightness of ~27 7 mag/D". This value 
is consistent at the \ .5cr level with the value of 26.7 7 mag/D" in- 
ferred by Zackrisson et al. (2009) for the resolved extragalactic 
background light emission. The isophotal radius of the stellar 
component of I Zw 1 8 and I Zw 1 8 C at //=26 B mag/D" was 
determined by P02 to be 8'.'8 and 5"5, respectively. From their 
combined isophotal area of ~340 D" and the above derived cps 
surface density of ~0.08 D"~ , we expect about 27 background 
cps in I Zw 18 and I Zw 18 C altogether. SBPs computed prior 



to and after removal of cps, were found in all bands to be con- 
sistent within lcr uncertainties, ensuring that cps contamination 
does not affect our conclusions in Sect. 3. 

SBPs were computed with the code iv (P02) (also referred 
to as Lazy by Noeske et al., 2006) that was specifically devel- 
oped for the study of irregular galaxies. This code permits a si- 
multaneous processing of co-aligned images of a galaxy in sev- 
eral bands and does not require a choice of a galaxy center, nei- 
ther does implicitly assume that the galaxy can be approximated 
by the superposition of axis-symmetric luminosity components. 
One of its key features is the computation of photon statistics 
within automatically generated irregular annuli that are adjusted 
to the galaxy morphology for each surface brightness interval 
H + Afi. This distinguishes code iv from other surface photome- 
try packages which generally employ ellipse-fitting to isophotes 
or photon statistics within elliptical annuli (e.g. meth. i of P96a, 
task ELLIPSE in IRAF, FIT/ELL3 in MIDAS), or approximate 
a galaxy by a single or several 2D axis-symmetric components 
(e.g. GIM2D, Simard (1998) and GALFIT, Peng et al. (2002)). 

3. Results 

3. 1. The photometric structure of I Zw 18 C 

In this section we investigate the photometric structure of 
I Zw 18 C, based on the combined HST ACS data in the filters 
V, R and 7. These allow to extend previous HST WFPC2 surface 
photometry (P02) by ^1 mag, with significantly reduced uncer- 
tainties below p. ~ 26 mag/D". In agreement with previous work, 
the SBPs of I Zw 18 C in all bands (Fig. 3a) were found to be 
nearly indistinguishable from one another, implying nearly con- 
stant radial colors. All SBPs display a narrow {R*<\") central 
excess and an outer (i"%R*%9") roughly exponential drop-off 
that mainly reflects the luminosity output from the host galaxy. 
A salient feature at intermediate radii (1"<R*S=3") is a shal- 
lower intensity increase than what inwards extrapolation of the 
outer profile slope to R*=0" predicts. An adequate fit to such 
SBPs therefore requires a modified exponential model that in- 
volves an extended (R*~3") central core of nearly constant sur- 
face brightness. One such fitting function (hereafter modexp), 
used by P96a to fit the host galaxy of BCDs (see Noeske et al., 
2003, for applications on near infrared (NIR) studies and a com- 
parison with the Sersic law) has the form: 



/(**) = /exp ■ [l - 61 eXpC-W))], 

where Pi(R*) is defined as 



R * \ 3 l R * i _ e 
P 3 (R*) = \ — + 



(1) 



(2) 



In addition to the central intensity 7q and scale length a of 
a pure exponential profile / exp = 7o • exp(R*/a), Eqs. 1&2 in- 
volve two further parameters: the central intensity depression 
e\ = A7/7o relative to an exponential model and the core ra- 
dius R c — 62 ■ a. The best-fitting modexp model to the V band 
SBP for ei,2=(2.4,0.8) (cf P02) is shown in Fig. 3a with the gray 
thick curve. It yields in all bands an extrapolated central sur- 
face brightness //e,o=21.7+0.2, a true central surface brightness 
po = Me.o - 2.5 log{\ - 62)=23.45 mag/D" and an a in the narrow 
range between 108 and 1 17 pc. The absolute V magnitude deter- 
mined from the modexp fit (-1 1.67 mag) corresponds to ~80% 
of the total luminosity of I Zw 18 C. 

Note that a direct determination of the intensity profile of the 
host galaxy of BCDs for radii R*< R c is generally prevented by 
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Fig. 3. (a) Surface brightness profiles (SBPs) of I Zw 18 C in V (F555W), R (F606W) and / (F814W). The thick gray curve shows a fit to the V 
SBP for R*>3'.'8 with the modified exponential fitting function Eq. 1 (modexp) for a core radius R c = 2.4 • a (dotted vertical line). The effective 
radius i? c ff and the radius Rgo enclosing 80% of the total V emission are indicated. The dashed line shows a linear fit to the outermost exponential 
part of the V SBP for R*>6", i.e. at the extremely faint (ji > 27 .6 mag/D") outskirts of the galaxy. The shadowed area corresponds to the lcr 
uncertainties of the V SBP. b) Radial V-R and R-I color profiles of I Zw 18 C, derived from the SBPs in the upper panel, c) Comparison of the 
best-fitting modexp model in panel a (thick gray curve) with the V SBPs of I Zw 18 C after partial removal of the brightest point sources (V*; 
squares). Open circles show the SBP of the unresolved stellar emission (Vd). computed after complete removal of compact (<075) sources with an 
unsharp masking technique, d) Comparison of the color profiles in panel b (open and filled circles) with the (V-I)a color profile (squares) of the 
unresolved stellar emission. The shadowed area depicts the lcr uncertainties of the V-R profile. 



the luminous young stellar component that typically dominates 
out to Rsp~2a (Papaderos et al., 1996b, hereafter P96b). Since 
young stellar clusters (SCs) can hardly be sufficiently resolved 
and subtracted out, even at the angular resolution of the HST, the 
chosen parameter set e;^ has to rely on plausibility arguments 
(see discussion in P96a and Noeske et al., 2003) and is to be 
considered approximative only. 

In the case of I Zw 18 C, however, due to the comparatively 
low surface density of young SCs, one stands a better chance to 
directly constrain the central form of the host's SBP. For this, 
we first subtracted out the brightest ~200 point sources from 
the galaxy using DAOPHOT (Stetson, 1979) and subsequently 
recomputed the V SBP from the residual emission. Figure 3c 
shows that this SBP (open squares, labeled V*) closely matches 
the best-fitting modexp model for intermediate to large R* . 
This agreement suggests that the adopted e\ o parameters yield 
a reasonable first-order approximation to the unresolved emis- 
sion of the host. SBP integration indicates that roughly 30% of 
I Zw 18 C's emission within R%q (Fig. 3) is due to point sources. 
This is rather a lower limit for the luminosity fraction of cps 
given that with the adopted procedure SCs could not be fully de- 



blended and subtracted out. This is apparent from the still strong 
(~2 mag) central peak of the V* SBP that is mainly due to in- 
complete removal of the central SC complex C and surrounding 
bright SCs (cf Fig. 4a). 

In an effort to better constrain the SBP of the unresolved 
stellar component of I Zw 18 C, we subsequently applied a flux- 
conserving unsharp masking technique (Papaderos et al., 1998, 
hereafter P98) to filter out all higher-surface brightness (HSB) 
clumpy features with <ff!5 and isolate the diffuse V emission 
only. The frames holding the compact (V c , R c and / c ) and diffuse 
(Vd, and /d) emission in each band are displayed in panels a 
and b of Fig. 4, respectively. The Vd SBP (open circles in Fig. 
3c) is at intermediate radii (l"^i?*^6") fairly comparable to the 
V* SBP, except for a nearly constant offset by «0.3 mag. Its 
corresponding absolute magnitude (-11.3 mag) is by about 0.75 
mag fainter than the integral value for I Zw 18 C (-12.05 mag) 
and can be regarded as characteristic for the unresolved stellar 
emission in the host galaxy. 

Note that all SBPs in Figs. 3a,c display at very faint lev- 
els (// ^ 27.6 mag/D") a shallower outer (R*^6") exponential 
slope (dashed line in Fig. 3a). This feature is certainly not due to 
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Fig. 4. (a) Three-color composite of the images R c , V c and / c (red, green and blue image channel, respectively), illustrating the spatial distribution 
of compact (<075) sources in I Zw 18 C. The regions C I, C II and C III defined by Izotov & Thuan (2004) are indicated. The blown up version 
of the central part of I Zw 18 C (lower- right) shows the ionized gas shell in the vicinity of the bright stellar cluster C (cf Dufour et al., 1996b). 
(b) Three-color composite of R4, Va and /j, of the unresolved stellar emission. The color coding is the same as in panel a. c) V-I color map of 
I Zw 18 C, revealing very blue colors (V-I <0) all over the southeastern third (region C I) of the galaxy. The contour corresponds to 25 V mag/D". 
d) Comparison of the spatial distribution of compact and point sources in V c and R c (blue and green channel, respectively) with the unresolved 
stellar background (/ d , red channel). 



point spread function (PSF) convolution effects since the maxi- 
mum extent of the ACS PSF at its lowest measured intensity (10 
mag below its central value) is ~1'.'5 (Jee et al., 2007). ne con- 
tamination from the main body can as well be excluded both on 
the basis of narrow-band imaging (e.g. Ostlin et al., 1996, P02) 
and because it would bluen V-I and R-I color profiles, in dis- 
agreement with the slight reddening of color profiles for R* ^6" 
(see below). This outermost SBP feature has evaded detection 
on previous surface photometry (P02) which, due to photometric 



uncertainties, was practically restricted to within the Holmberg 
radius. Because of its faintness (merely 3% of the total emission) 
and low surface brightness (~27.6 - 29 mag/D"), its reality can 
not be established beyond doubt even with the present data. If 
due to an underlying stellar host of constant exponential slope to 
R*=Q", its no 23.6 mag/D") and a (160 pc) would qualify it 
as an LSB dwarf with a My ^-l 1 mag (equivalent to ~38% of 
I Zw 18 C's emission). 
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We turn next to the color distribution in I Zw 18 C. Figure 3b 
shows that for 1"£R*£6" both the V-R and R-I index is very 
blue and nearly constant (=0 with a standard deviation about the 
mean of 0.05 mag). Since R*-6" encompasses practically the 
total emission of I Zw 18C, this color may be regarded as rep- 
resentative for the galaxy as a whole. In the outermost periph- 
ery of the galaxy, where the shallower outer exponential compo- 
nent appears, color profiles hold a hint for a slightly redder R-I 
(0.16+0.09) without a notable change in V-R. The color profiles 
in their innermost (R*<1") part are dominated by luminous cps 
in the surroundings of region C (cf Fig. 1) and show a large scat- 
ter (up to 0.15 mag) around mean values of and 0.06 mag for 
V-R and R-I, respectively. 

In order to place meaningful constraints on the evolution- 
ary status of I Zw 18C, it is necessary to verify that these blue 
colors are not due to the luminosity-weighted average of bright 
blue SCs with a red underlying stellar background. The latter 
could readily escape detection, even when dominating the stellar 
mass (see e.g. P98). This concern is underscored by two atypical 
properties of I Zw 18 C: First, the galaxy exhibits a weak color 
gradient along its major axis (cf Aloisi et al., 1999; Izotov et al., 
200 1 a), with its northwestern tip being redder (B - V = 0.05 mag, 
V-I = 0.2 mag) than the southeastern one (B-V = -0.07 mag and 
V - I = -0.2 mag). This is apparent from Fig. 4c, from which a 
V-I color as blue as -0.2 mag can be read off all over the south- 
eastern third of the galaxy (region C I in the denomination by 
IT04). By contrast, the average V-I color in the central (C II) and 
northwestern (C III) part of the galaxy is redder (>0 mag) with 
several local color maxima (>0.3 mag) associated with cps. It 
is thus conceivable that the bluer southeastern and redder north- 
western galaxy half counter-balance each other, thereby intro- 
ducing an overall blue mean radial color. Another characteristic 
of I Zw 1 8 C that could additionally conspire in diminishing ra- 
dial color gradients is that the surface density of its cps tends 
to be spatially anti-correlated with the unresolved stellar back- 
ground of the host. This is illustrated in Fig. 4d where the diffuse 
emission of the latter (Id', red channel) is overlaid with the cps in 
the V c and R c images (blue and green channel, respectively): it 
can be seen that the diffuse component peaks at the northwestern 
half of I Zw 18 C (C II and C III) in which it accounts for =50% 
of the V line-of-sight intensity, whereas its contribution drops to 
$30% in the southeastern half of I Zw 18 C (region C I) where 
most of the bright blue cps are located. Consequently, especially 
in region C I, a hypothetical red stellar background could readily 
escape detection on radial color profiles. 

From such considerations, and in order to infer the colors 
of the unresolved stellar component of I Zw 18 C in an as un- 
biased manner as possible, we additionally computed color pro- 
files based on Vd, Rd and Id SBPs only. In all cases, we found a 
good agreement with the results initially obtained from the total 
emission with the exception of a central (R*H") red peak with 
mean values of (V-R) A =Q- 1±0.03 and (V-/) d =0. 16+0.05. This 
color excess might be attributed to a stellar age gradient or en- 
hanced extinction in region C and surroundings. Note that Izotov 
et al. (2001a) estimated from spectral synthesis models an ex- 
tinction coefficient C(Hy6)=0. 1-0.3 for region C, corresponding 
toA^O.2-0.65 mag. 

As apparent from panel d of Fig. 3, the V-I profile of the un- 
resolved stellar component ((V-I)d, squares) is fairly compara- 
ble to the V-R and R-I profiles (panel b), revealing a nearly con- 
stant color of =0+0.04 mag within 1 " <R* <6" and a slightly red- 
der value (0.2+0.08 mag) in the extreme periphery of the galaxy. 
We are therefore led to conclude that the overall blue colors of 
I Zw 18 C for l"^jR*<6"are not dictated by bright young SCs 
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Fig. 5. (upper panel) V, R and / SBPs of I Zw 18. The thick-solid 
and thin-dashed lines show linear fits to the V SBP in the radius range 
7"<R*<15" and R*>7", respectively. The effective radius R e ff and the 
radius R 80 enclosing 80% of the total V emission are indicated. The V 
SBP derived by P02 after subtraction of the HyS and [Om],U4959,5007 
emission lines (referred to as V"), shifted by +1 mag for the sake of 
clarity, is included for comparison, (lower panel) V-R (filled circles), 
V-I (squares) and R-I (open circles) color profiles, computed from the 
SBPs in the upper panel. The vertical bar indicates the V-I color range 
inferred by P02 from ground based and HST PC data for region at at 
the southeastern tip of I Zw 18 (V-I(u)) = 0.17 ...0.32 mag). Since 
nebular emission is relatively weak in that region, its color can be con- 
sidered representative for the stellar host of I Zw 18. The mean V-R 
color of I Zw 18's host after subtraction of nebular line emission (V- 
R' =0.12+0.04 mag; P02) is indicated by the horizontal line. Note that 
the colors V-I(<±>) and V'-R' of the stellar host galaxy are by, respec- 
tively, S0.8 mag redder and »0.4 mag bluer than the colors of the expo- 
nential nebular envelope (R*>6"). 



spread all over the body of the galaxy but as well characteristic 
for its unresolved stellar component. 

3.2. The photometric structure of I Zw 18 

The V, R and / SBPs of I Zw 18 (Fig. 5) compare well to those 
presented by KO00 and P02. They exhibit a central (R*<6") 
high-surface brightness core and a nearly exponential LSB en- 
velope extending out to R*~2Q". Profile fitting in the radius 
range 7"<R*< 15" (solid line) yields for the LSB component 
an a = 2'.'3 (210 pc) and a luminosity fraction of =30%. At 
fainter levels 27.6 - 30 mag/n") all SBPs exhibit a shallower 
exponential slope with a = 471 (370 pc). This SBP feature 
is attributable to the faint filamentary emission in the extreme 
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Fig. 6. V-R (left panel) and R-I (right panel) color map of I Zw 18, displayed in the range between and 0.4 mag and -0.9 and mag, respectively. 
Crosses mark the star-forming regions NW and SE. The region termed a> by P02 (cf Fig. 1) is indicated at the southeastern part of the image. The 
insets in the left panel show a magnified view of a shell «5" northwards of region NW. 



northwestern and southeastern periphery of I Zw 18 (cf. Fig. 8) 
that has evaded detection on previous HST WFPC2 imagery and 
which accounts for no more than 1 % of the total luminosity. A 
linear fit to the whole LSB envelope (7"<R*<20"; dashed line) 
yields a mean a = 2'.'9 + 07 13 (-270 pc). 

The radial color profiles of I Zw 18 (lower panel of Fig. 5), 
reflect the increasing contribution of ne to the observed line-of- 
sight intensity with increasing galactocentric radius, in agree- 
ment with previous evidence (P02). The V-R color (filled cir- 
cles) increases roughly linearly from S50 at R*=Q" to 0.55 mag at 
R*=9" where it levels off to a nearly constant value. The corre- 
sponding, relatively strong color gradient of y + =0.6 mag kpc -1 
is rather typical for evolved BCDs (P96b), suggesting, at first 
glance, that SF activity in I Zw 18 occurs within a more ex- 
tended, old underlying host. This interpretation is, however, im- 
mediately challenged upon inspection of the extremely blue V-I 
and R-I profiles. The latter display already in their inner part 
(2"5=R*<6") mean values of -0.21 mag and -0.44 mag, respec- 
tively, both inconsistent with the red V-R color, if the emission is 
assumed to be dominated by stars. More impressively, at R*~6", 
i.e. roughly at the transition radius between the HSB core and the 
exponential LSB envelope of SBPs, either color index shows a 
sudden decrease to values as blue as V-/=-0.61+0.13 mag and 
R-/=-l.l+0.08 mag, which then remain nearly constant out to 
R*~2Q". 

That such colors can not be of stellar origin is apparent al- 
ready from the fact that, even for an 05 V star, the V-I and R-I 
color is -0.32 mag and -0.18 mag, respectively, i.e. 0.3 to 0.9 
mag redder than the color of the LSB envelope of I Zw 18. More 
generally, as pointed out by P02, there is no stellar population, 
regardless of star formation history, age and metallicity that can 
reproduce the observed combination of red (~0.6 mag) V-R with 
blue (-0.6 ... -1 . 1 mag) V-I and R-I colors of the LSB envelope 




Fig. 7. V-I map of the central 15'.'6xl5'.'6 of I Zw 18, displayed in the 
color range between -0.5 and 0.15 mag. The contours are computed 
from the R' image by P02 (see also Fig. 8) and go from 19 to 25 mag/D" 
in steps of 0.5 mag. Note the very blue (<-0.5 mag) rim encompass- 
ing the northwestern star-forming region. The mean V-I color over the 
southeastern region a> is determined to be 0.23 mag with an upper value 
of 0.3 mag in its reddest quarter. 

of I Zw 18. Quite to the contrary, such colors can solely be ac- 
counted for by ne. 

Indeed, comparison of Fig. 5 with Fig. 12 of P02 reveals that 
the radius R*=6"- 7" where the steep V-I and R-I color drop- 
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Fig. 8. Composite image of I Zw 1 8 and I Zw 1 8 C with the R, I and V bands shown in the red, green and blue channel, respectively. The image 
discloses a complex network of ionized gas filaments extending as far out as 2.6 kpc, approximately twice the distance reported from previous 
HST WFPC2 studies, corresponding to ~16 exponential scale lengths a of the stellar host galaxy of I Zw 18 (160 pc; P02). The contours, adapted 
from P02, are computed from HST WFPC2 R data, after two-dimensional subtraction of the Ho- emission (referred to as R'), they thus delineate 
the morphology of the stellar and nebular continuum emission in the main body of I Zw 18. Contours go from 19 to 25 R' in steps of 0.5 mag. 
Note that I Zw 18 (i.e. its stellar component, as depicted by the R' contours) and I Zw 18 C show a remarkably similar structure. 



off occurs is identifiable with the radius where the line-of-sight 
contribution of stellar emission steeply decreases from a plateau 
value of -40% for 3£R*£6" to less than 20%. Therefore, the 
much deeper data studied here confirm and strengthen the previ- 
ous conclusion that the extended (6"^R* ^20") exponential LSB 
envelope of I Zw 18 is due to ne. 

The gaseous nature of the LSB envelope of I Zw 18 is also 
evident from the color maps in Figs. 6 and 7. These show a clear 
correspondence to radial color profiles, most notably a strong 
core-envelope color contrast (0.5 - 1.5 mag) and remarkably 
uniform colors for the envelope (V—R~ 0.5 ...0.6, R—Isn—1A 
mag, V-I~-l mag) over a spatial scale of 9 kpc 2 . The northwest- 
ern super-shell, for example, though prominent on direct images 



(cf Fig. 1), is barely distinguishable from its gaseous surround- 
ings on V-I and R—I color maps, suggesting a nearly constant 
spectral energy distribution (SED) all over the LSB envelope. 

The relative extent of the stellar component that is confined 
to the compact HSB core with respect to the nebular LSB enve- 
lope is better illustrated in Fig. 8. The image reveals a complex 
network of ionized gas filaments extending as far out as 2.6 kpc, 
twice the galactocentric distance previously reported from HST 
WFPC2 studies and equivalent to ~16 exponential scale lengths 
a of the stellar host of I Zw 18 (160 pc; P02). The contours, 
adapted from P02, are computed from HST WFPC2 R data, af- 
ter two-dimensional subtraction of Ha line emission (referred 
to as R'), they thus delineate the morphology of the stellar and 
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Fig. 9. Three-color image of I Zw 18 displaying the unresolved stellar emission in the red channel and the R and V emission of compact 
ne-dominated regions in the green and blue channel. White areas in the central (red) part of the galaxy depict regions whose colors are strongly 
affected by nebular emission, as apparent from their extremely blue V-I and R-I (< -0.4 . . . -0.8) colors. The inset, showing a magnified portion of 
the western super-shell, is meant to illustrate a subtle spatial shift between the V c and R c emission. This can be attributed to a spatial displacement 
by 30 - 100 pc between the intensity maxima of the [Olll]/1/I4959,5007 and Ho- emission lines. 



nebular continuum emission. Note that I Zw 18 (i.e. its stellar 
component as depicted by the R' contours) and I Zw 18 C show 
a remarkably similar structure. 

3.2.1 . The impact of nebular emission on the colors of the 
stellar component of I Zw 1 8 

We discuss next in more detail the effect of ne contamination 
in the HSB core, i.e. the region subtended by the 25 R' mag/D" 
isophote in Fig. 8. Our color maps reveal here a substantial sub- 
structure that was not sufficiently recovered on previous data. A 
previously known feature is the ne rim (EW(Hor)= 1500-2000 A; 
Ostlin et al., 1996; Vflchez & Iglesias-Paramo, 1998; Papaderos 



et al., 2001; Izotov et al., 2001a) that encompasses the SF region 
NW. Its V-I and R-I colors are much bluer (-0.8 mag and -0.6 
mag, respectively) than those of the centrally located young SCs 
(-0.15 mag and -0.23 mag, measured within a 2'.'8x2'.'8 aper- 
ture) where the EW(Ha) is much lower (~200 A; Izotov et al., 
2001a). 

This striking spatial anti-correlation between stellar surface 
density and emission-line EWs, with the V-I color steeply de- 
creasing with the increasing EW in the periphery of ionizing 
SCs is essentially identical to that described in the XBCD SBS 
0335-052E by P98: the bluest V-I colors (-0.8 mag) in that 
galaxy were not observed at the position of its young SCs but 
in their periphery, over an extended horseshoe-shaped gaseous 
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rim some 500 pc offset from the former. Several similar exam- 
ples are documented among XBCDs and BCDs (e.g. Papaderos 
et al., 1999; Guseva et al., 2001; Fricke et al., 2001; Ostlin et al., 
2003; Guseva et al., 2004; Papaderos et al., 2008) both on small 
and large spatial scales. 

Of considerable interest is also the small-scale contamina- 
tion of colors in the HSB core of I Zw 18 due to ne. To better il- 
lustrate its impact, we computed a pseudo three-color composite 
(Fig. 9) using a combination of unsharp masked images (see dis- 
cussion in Sect. 3.1). In order to partly suppress the stellar con- 
tribution in V c and R c images (blue and green channel, respec- 
tively), we subtracted from the latter a pseudo stellar continuum 
using a scaled version of the I c image. This procedure allowed 
to better isolate and visualize regions where ne contamination is 
strongest. The red channel of Fig. 9 holds the unresolved emis- 
sion in the / band (/j), it thus primarily reflects the stellar host 
galaxy. In the resulting three-color overlay, white patches on top 
the reddish background of the HSB core depict regions where ne 
dictates colors. Just like in radial profiles and color maps (Figs. 
5, 6 and 7), the footprint of ne is a combination of red V-R with 
extremely blue V-I and R-I colors. Some examples include the 
regions labeled 1 through 8, whose respective colors are (0.1, 
-0.35, -0.45), (0.17, -0.58, -0.76), (0.05, -0.41, -0.45), (0.23, 
-0.43, -0.67), (0.24, -0.45, -0.69), (0.21, -0.53, -0.74), (0.12, 
-0.35, -0.47) and (0.23, -0.23, -0.46). The total area of these 
severely ne-contaminated regions of «30 □" is equivalent to the 
isophotal size of I Zw 18 at 22 mag/d" (32 □") and to 20% of 
that of its stellar component (the region subtended by the 25 R' 
mag/D" isophote in Fig. 8). 

It is noteworthy that ne shows significant substructure on 
spatial scales of a few pixels, with little spatial correlation to the 
local stellar background. Therefore, its treatment and subtrac- 
tion as uniform foreground emitting layer in CMD studies could 
result into systematic uncertainties that might not be fully ac- 
counted for by the standard error budget. Subtle spatial displace- 
ments between the intensity maxima of strong nebular emission 
lines may be of some importance in this regard. An example is 
given in the inset of Fig. 9 where we show magnified a portion 
of the western super-shell. It can be seen that the inner layer 
along the shell appears blueish whereas the outer one green- 
ish. This can be plausibly attributed to a displacement by 07 1 
- 073 between the [Om]/U4959,5007 and Ha emission lines 
which are registered in V (blue channel) and R (green channel), 
respectively. Another illustrative example is given for a shell 
~5" northwards of region NW (Fig. 6) whose inner (southeast- 
ern) and outer (northwestern) layer differ by ^0.3 mag in their 
V-R color, pointing to a variation by ^30% in the [Ollll/Ha 
ratio across the shell. Note that systematic variations of the 
[O Ill]/15007/Hy6 ratio by a factor ~3 on spatial scales of a few 
ten pc have been documented by integral field unit spectroscopy 
in a number of Hll regions and BCDs (see e.g. Westmoquette et 
al., 2007; Cairos et al., 2009; Relano et al., 2010; Monreal et al., 
2010). 

Whereas spatial displacements between nebular lines have 
certainly no measurable effect on CMD studies of galaxies with 
faint ne, they may be of some relevance in the case of I Zw 18. 
This is because the local background level in point source pho- 
tometry studies may be overestimated by a different amount in 
different bands, depending on the position of a star relative to 
a shell and the specifics of the local background determination. 
Figure 10 illustrates schematically two special cases: for a star 
located near the inner interface of the shell (case a) the circular 
annulus within which the local background is determined cap- 
tures the inner [Olll] -enhanced portion of the shell and misses its 




100 pc 



Fig. 10. Schematic illustration of the effect that the slight spatial dis- 
placement between the [Om]/i/i4959,5007 and the Ha emission lines 
across an expanding shell of ionized gas might have on V and R photom- 
etry of stars in its close vicinity. The shell, of thickness 100 pc (=al"at 
the distance to I Zw 18) comprises an inner and outer zone with, respec- 
tively, a higher and lower [OlIl]/Ho' ratio (cf the discussion related to 
Figs. 6,left and 9). Two cases may be distinguished: in the case of a star 
located close to the inner boundary of the shell (a), the local V back- 
ground within a circular annulus will be overestimated due to the en- 
hanced contribution of the [Olll] lines in the inner shell interface. This 
would result in an artificially red V-R color. The opposite would be the 
case for (b) where the local R background would be over- subtracted due 
to the lower [Om]/Ha ratio in the outer layer of the shell, leading to an 
artificially blue V-R color. The principal effect would therefore be that 
the magnitudes of stars in the vicinity of shells would be underestimated 
by a different amount in different bands, with the V-R colors of stars in 
the interior of the shell being redder than in its outer part. 



outer Hcr-enhanced layer. This could lead to an over-subtraction 
of the local V background, hence, to a reddening of the V-R 
color. The opposite might be the case for a star close to the outer 
boundary of the shell (case b), whose color would appear bluer 
due to the over-estimation of the R (Ha enhanced) background. 
In the simple geometry of Fig. 10, the principal effect would then 
be a redder color for stars in the shell interior and vice versa, in 
addition to underestimated stellar magnitudes because of over- 
subtraction of the local background. An examination of the cu- 
mulative effect that multiple overlapping ne shells could have on 
CMD studies of I Zw 18 might therefore be of some interest. 

Of special relevance to the study of the evolutionary sta- 
tus of I Zw 18 (Sect. 4.2) are the colors of region at (Fig. 
1). This region, roughly delimited by the rectangular area at 
(RA,DEC)=(578 . . . 476,-676 . . . -474) relative to NW (see Figs. 
6 and 7) is relatively free of ne contamination (cf the EW(Ho') 
map in Fig. 7 of Izotov et al., 2001a), its colors can therefore be 
used to place constraints on the age of I Zw 18. From the present 
data we infer a mean V-R color of -0.06 mag and a R-I and V-I 
color of xQ.2 mag with values of V-R - R-I ~ 0.15 mag and 
V-I=0.3 mag within its reddest quarter. The V-I color range in- 
ferred for region a> from HST ACS data is in good agreement 
with the values 0. 17 - 0.32 mag previously determined by P02. 



12 



I Zw 18 as morphological paradigm for rapidly assembling high-z galaxies 



4. Discussion 

4. 1 . Constraints on the evolutionary status of I Zw 18 C 

In Fig. 1 1 we compare the observed colors of I Zw 18 C with the 
predicted color evolution for a stellar population forming instan- 
taneously (SFH1, dotted line), with an exponentially decreasing 
star formation rate (SFR) and an e-folding time t — 1 Gyr (SFH2, 
solid line) and continuously with a constant SFR (SFH3, dashed 
line). The theoretical curves were computed with the evolution- 
ary synthesis code Pegase 2.0 (Fioc & Rocca-Volmerange, 1997) 
for constant stellar metallicity of Z=0.001 and Salpeter initial 
mass function (IMF) between 0.1 and 100 M Q , and do not in- 
clude nebular emission. 

As already pointed out in Sect. 3.1, the blue and nearly con- 
stant (0+0.05 mag) V-I and R-I colors of I Zw 18 C down to 
27.6 mag/D" imply that the photometrically dominant stellar 
population in this system is almost uniformly young. Formal up- 
per age estimates within ler uncertainties range from t\ — 15 
Myr for SFH1 to t 2 -35 Myr for SFH2.3. The latter SF scenar- 
ios are, however, incompatible with the data since star formation 
starting at t 2 and continuing to the present would have given rise 
to amble ne with an EW(Hor)^800 A. The uniformly blue col- 
ors of I Zw 18 C, in connection with the absence of ne every- 
where but in its central region indicate that SF activities must 
have ceased very recently, some ~20 Myr (= ? c ) ago. If so, the 
estimated SF duration (h — t c - t\ ) would imply an almost in- 
stantaneous SF process, rapidly synchronized over the projected 
area of the galaxy (~1 kpc 2 ) at probably supersonic speeds. 

A conceivable scenario might invoke sequential star forma- 
tion from the redder northwestern tip towards the bluer south- 
eastern tip of I Zw 18C(see also Aloisi et al., 1999; Izotovet al., 
2001a). The respective colors of those regions translate by SFH1 
to an age difference of tsf =40-80 Myr. Taking tsf ~ 60 Myr 
as an indicative time scale for the spatial progression of SF ac- 
tivities along the projected major axis of I Zw 18 C (1 kpc), one 
can estimate an average SF propagation velocity of msf ~ 20 km 
s _1 . The latter is comparable to the msf inferred for the XBCD 
SBS 0335-052E (~ 20 - 35 km s" 1 , P98, Reines et al., 2008), of 
the order of the sound speed in the warm ISM. Induced gas col- 
lapse along the collisional interface of a super-shell expanding 
from the northwestern tip of I Zw 18 C between 80 Myr (rsp+t c ) 
and t c ago might offer a tenable, though not necessarily unique 
explanation for propagating star formation: following McCray 
& Kafatos (1987), the radius (pc) of a SF-driven super-shell 
can be approximated as 

R sh = 269 (^)' /5 tf (3) 

where L m ^ is the mechanical luminosity injected into the ISM 
by stellar winds and SNe in 10 38 erg s" 1 , «o me ambient gas den- 
sity in crrT 3 and tq the dynamical expansion time in 10 7 yr. A 
rough estimate on the mean mechanical luminosity can be in- 
ferred from the absolute V magnitude of I Zw 18 C (-12.2 mag) 
which translates for continuous star formation over tsf and an 
ensuing quiescent phase over the past t c to a stellar mass of 
w 2.1 x 10 5 M Q and a mean SFR of 3.5x 10~ 3 M Q yr 1 . The mean 
£m,38 over t sf may be estimated from Starburst99 (Leitherer et 
al., 1999) to be 14.6. Equation 3 then yields for R sh = 1 kpc 
and t-j — 6 (= tsf), an ambient gas density no * 4 cm -3 . This 
value does not seem to be unrealistic, given that the VLA HI 
map by van Zee et al. (1998) indicates for I Zw 18 C a Hi sur- 
face density of ^ 4 x 10 20 cirT 2 , which after inclusion of the 
helium contribution translates for a disk thickness of ~50 pc to 
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Fig. 11. Comparison of the V-I and R-I colors of I Zw 18 C and of re- 
gion to in I Zw 1 8 with model predictions for a stellar population form- 
ing instantaneously (SFH1: thick dotted curve) or continuously, with 
an exponentially decreasing or constant star formation rate (SFH2 and 
SFH3: thick solid and dashed line, respectively). The models have been 
computed with Pegase 2.0 (Fioc & Rocca-Volmerange, 1997) for a con- 
stant metallicity (Z=0.001) and a Salpeter initial mass function between 
0.1 and 100 M e , and do not include nebular emission. Thin curves cor- 
respond to the same star formation histories but assume an IMF trun- 
cated above 5 M Q . The vertical bars labeled IZw 18to> depict the range 
between the mean and reddest color in region u> of I Zw 18 (V-/=0.2 
. . . 0.3 and fi-/=-0.06 . . . 0.15). Bars labeled I Zw 18 C (l"-6") indicate 
the mean color of I Zw 18 C in the radius range 1"<R*<6" («0±0.05 
mag). The mean colors and their lcr uncertainties in the extreme pe- 
riphery (R*>6") of I Zw 18 C at ji i 27.6 mag/D" are depicted by the 
rectangular areas. 



«o ~ 3 crrr 3 . As the passage of the SF front along I Zw 18 C has 
likely been accompanied by the depletion of its molecular con- 
tent and the dispersal of its ISM, the present gas surface density 
should rather represent a lower estimate on the average «o over 
Tsf- Clearly, the considerations above rely on strongly simpli- 
fying assumptions (constant L m ^% and no, no radiative losses in 
the shell, coplanar face-on geometry) and are meant to be indica- 
tive only. They merely touch upon one among several possible 
scenarios behind propagating star formation in I Zw 18 C (e.g., 
gravitational interaction with the main body or gas infall from 
the Hi halo) and are invoked to demonstrate that this process can 
account for some important characteristics of the galaxy. These 
include its uniformly blue colors with merely a weak color gradi- 
ent along its major axis, the virtual absence of ne and the higher 
surface density of the unresolved stellar background in its north- 
western half as possible signature of the disintegration of its old- 
est SCs. However, propagating star formation between tsf+^c 
and t c can alone not explain the slightly redder color of stellar 
populations in the outskirts of I Zw 18 C. 
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4.1.1. Stellar diffusion and its effect on age 

determinations in the faint periphery of I Zw 18 C 

In the following we further explore the evolutionary history of 
I Zw 18 C by considering the properties of its faint (27.6 - 
29 mag/rj") stellar periphery. The (V-I) d (0.2+0.08) and R-I 
(0.16+0.09) color of the latter translate by SFH2 to an age of 
~130 Myr with an upper bound of ~500 Myr at the lcr level, 
if purely stellar emission is assumed. For models including ne, 
age estimates would rise to 270-900 Myr, depending on the 
color considered. However, even for those high ages, this SFH 
model (and SFH3 alike) predicts an EW(Ha) between 150 A 
and ^300 A in clear conflict with the absence of ne in the out- 
skirts of I Zw 18 C. Models invoking continuous star forma- 
tion to the present are thus fundamentally incompatible to the 
data. Cessation of SF activities over the past t c would alleviate 
the contradiction, it would, however, at the same time imply a 
steeper color evolution, thus a younger age, in better agreement 
with the lcr upper estimate ?sfhi ~ 100 Myr read off Fig. 1 1 for 
the instantaneous star formation model. 

As we shall argue next, stellar diffusion and the resulting ra- 
dial stellar mass filtering effect described by P02 could have a 
significant effect on the color distribution of young galaxy can- 
didates and add an important element towards a consistent evolu- 
tionary picture for I Zw 18 C. Already at a mean radial velocity 
of w r ^ 4 km s _1 , less than the velocity dispersion of the neu- 
tral gas (o" H i ~6-8 km s~'), a star born in the central part of 
I Zw 18 C could migrate within rdiff = tsf + t c (80 Myr) out 
to r^ff ~ 300 pc from its initial locus. This galactocentric dis- 
tance corresponds to the semi-minor axis of I Zw 18 C at 29 
mag/D", i.e. it is topologically identifiable with the redder, outer 
exponential SBP feature in Fig. 3. A consequence of stellar dif- 
fusion with the above radial velocity is that any stellar genera- 
tion reaching out to r&s can not be bluer than 0.3 mag in V-I 
(approximately the lcr color bound, depicted by the rectangu- 
lar area in the upper panel of Fig. 11), since on its way to that 
radius it will have been de-populated from stars with lifetimes 
^ fdiff (or, correspondingly, from stars more massive than Mdiff). 
The commonly adopted SFH2.3 models would then inevitably re- 
sult in an overestimation of the stellar age in the LSB periphery 
(R*-rdift) of young galaxy candidates, such as I Zw 18 C. This 
is because such SFH parametrizations are throughout applied on 
the implicit assumption that stellar populations form in the ra- 
dial zone where they are currently observed, hence they invari- 
ably contain a certain fraction of massive blue stars younger than 
Tdiff. Obviously, a match between predicted and observed colors 
is then only possible for an age fsFH2,3 exceeding the true stellar 
age f* at the radius R*, i.e. the condition 

Tdiff ^ fSFHl < f* < fSFH2,3 (4) 

applies throughout, in particular when models including ne are 
adopted. 

That, in the presence of stellar diffusion, the age of the LSB 
host of young BCD candidates may be overestimated when con- 
tinuous star formation models are used is especially apparent 
when projection effects are additionally taken into account: in a 
spherical-symmetric geometry, the colors registered at R*= rdiff 
reflect the luminosity-weighted average of stars at radii > rdiff, 
since the observer's line of sight crosses more remote foreground 
and background galaxy zones of longer Tdiff. Therefore, even 
though SFH2.3 may provide a reasonable approximation to the 
integral SFH of some star-forming dwarf galaxies (SFDGs), it 
should not be taken for granted that they are universally applica- 
ble to spatially resolved age-dating studies. This is already sug- 



gested by the failure of e.g. SFH3 to reproduce the generally red 
colors (V-7=0.9, B-R=l.l) of the host galaxy of typical BCDs 
(see e.g. P02), even for ages exceeding the Hubble time. 

The effect of stellar diffusion to the radius rdiff is equivalent 
to in situ star formation at rdiff with an IMF truncated at masses 
above a radially decreasing limit Mdiff (rdiff) (hereafter tIMF). To 
illustrate this point, and as a minimum consistency check, one 
may consider continuous star formation according to SFH2,3 with 
a tIMF limited to stars with a lifetime ^ Tdiff. This timescale 
roughly corresponds to the main sequence lifetime of a B4 star 
(94 Myr), placing an upper cutoff of 5 M to the tIMF. The color 
evolution for SFH2.3 with that tIMF is depicted in Fig. 1 1 with 
thin lines. It can be seen that SFH2,3+tIMF models reproduce 
the observed colors at an age of ~130 Myr with an lcr upper 
bound of 150-250 Myr and consistently account for the absence 
of nebular emission. 

A conceivable alternative to stellar diffusion is in situ star 
formation in the LSB periphery of I Zw 18 C in conjunction with 
stochastic effects on the IMF (see Cervino & Mas-Hesse, 1994, 
and references therein). The latter are to be expected when low 
masses are involved in the SF process, as is typically the case for 
SFDGs, and to affect the sampling of the IMF in its intermediate- 
to-high mass range. Cervino & Mas-Hesse (see their Fig. 1) ar- 
gue, however, that stochastic effects do not result in a truncated 
IMF above a certain stellar mass, but some massive stars always 
form. On the statistical average, one would therefore expect neb- 
ular emission to be present in the outskirts of I Zw 18 C, in dis- 
agreement with the observations. It should be noted, though, that 
this conclusion is likely dependent on the detailed physics of star 
formation, from which the high-mass end of the IMF is popu- 
lated and the time scales for low and high mass star formation. 
With certain assumptions (see e.g. Weidner & Kroupa, 2005; 
Pflamm-Altenburg et al., 2009) a truncated IMF can be real- 
ized in peripheral galaxy zones, lifting the discrepancy described 
above. There might therefore be variants of the traditional static 
scenarios, i.e. scenarios implicitly assuming that stars are born in 
the radial zone where they are observed, that can quantitatively 
reproduce the radial color distribution of I Zw 18 C and SFDGs 
in general. It is not clear, however, whether such static scenarios, 
and their assumptions on the radial dependence of the IMF, can 
consistently account for the specific structural characteristics of 
SFDGs, and in particular BCDs, such as e.g. the exponential- 
ity of the stellar LSB host and the confinement of star-forming 
activities to within R*^ 2a (see e.g. P96b). 

Stellar diffusion and its role on the buildup of SFDGs consti- 
tutes an unexplored territory in dwarf galaxy research. A caveat 
of the above discussion is that it assumes young stars to initially 
drift apart freely at a roughly constant velocity u r (^ cr HI ) over 
~ 10 8 yr. Clearly, the validity of this simplifying assumption 
needs to be investigated, and the initial kinematics of the newly 
formed stars as well as the form of the gravitational potential be 
included in the analysis. High-velocity resolution (S 10 km s~') 
integral field spectroscopy might provide key insights into the 
velocity patterns and the possible outwards migration of bound 
and unbound stellar clusters in young starbursts. A closer in- 
vestigation of the effects of diffusion would be of considerable 
interest, among others because it offers a mechanism that natu- 
rally produces radial color gradients in galaxies. As such, it may 
be of relevance for a range of topics, such as e.g. the Red Halo 
phenomenon (Bergvall & Ostlin, 2002; Zackrisson et al., 2006; 
Bergvall et al., 2010). 

In summary, our results suggest that the overall photomet- 
ric properties of I Zw 18 C can be consistently accounted for 
by a single episode involving solely sequential star formation 
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(wsf ~ 20 km s -1 ) and stellar diffusion (u r « 4 km s -1 ) over the 
past ~ 100-200 Myr. Whereas an ultra-faint substrate of ancient 
stars can not per se be ruled out, the lines of evidence presented 
here indicate that I Zw 18 C is a cosmologically young object. 
This conclusion concurs with that by Izotov et al. (2001a) who 
inferred an upper age of ~100 Myr for I Zw 18 C from spectral 
synthesis models. Jamet et al. (2010) reached a similar conclu- 
sion, that I Zw 18 C is younger than ^125 Myr, from a proba- 
bilistic analysis of HST ACS CMDs. 

4.2. Constraints on the evolutionary status of I Zw 18 

We turn next to the evolutionary status of I Zw 18's main body. 
The V—I and R-I color range in region u> (Fig. 11), translate 
by SFH2.3 to an age between ~100 and $250 Myr. The possi- 
bility of I Zw 1 8 and I Zw 1 8 C forming a co-evolving pair of 
dwarf galaxies that underwent a roughly synchronous strong re- 
cent evolution can not be ruled out. 

By considering the reddest quartile of region o>, one obtains 
from SFH2.3 an upper age of ~500 Myr. The age span inferred 
from the present study is consistent with that obtained from the 
mean B'-V (0.09+0.04 mag) and V'-R' (0.12+0.04 mag) for 
the stellar host galaxy of I Zw 18 after two-dimensional sub- 
traction of strong nebular lines from broad band images (P02). 
Note, however, that the latter emission-line free colors are likely 
slightly overestimated, as they do not include corrections for the 
red (fi-V=0.34, V-R=0.64, see e.g. Kriiger et al., 1995) nebular 
continuum. Notwithstanding this fact, even when taken within 
their +lcr bounds at face value, they imply by SFH2 an age of 
~0.8 Gyr, translating into a mass fraction of ^50% for stars 
younger than 0.5 Gyr. Evidently, this conclusion (see also Hunt 
et al., 2003) is not in conflict with the presence of a small num- 
ber ($20) of stars with ages between 0.5 and ~1 Gyr in CMDs 
(Aloisi et al., 1999; Ostlin, 2000; Ostlin & Mouhcine, 2005). On 
the other hand, it is important to bear in mind that estimates on 
the mass fraction of young stars depend critically on the adopted 
SFH. 

4.2.1 . An extended underlying stellar disc in I Zw 1 8? 

We next revisit the hypothesis envisaged by Legrand (2000), ac- 
cording to which the formation of I Zw 1 8 is occurring through- 
out its Hi halo (60"x45") at an extremely low SFR over the past 
~14 Gyr. This model can reconcile the low gas-phase metallicity 
of I Zw 18 with a cosmological age, and predicts a high degree of 
chemical homogeneity in its warm ISM. From the photometric 
point of view, the main prediction from the Legrand model is an 
extended ultra-LSB underlying stellar disk (Ji ^28 V mag/rj"). 
In a subsequent study, Legrand et al. (2001) have generalized the 
scenario of 'slowly cooking' dwarfs to BCDs, arguing that their 
main metallicity and stellar mass contributor is continuous star 
formation, rather than coeval starbursts. 

The hypothesis of a stellar disk beneath the nebular enve- 
lope of I Zw 18 was later on investigated by P02 on the basis 
of emission-line free HST WFPC2 B', V' and R' images. By 
consideration of photometric uncertainties, P02 argued that the 
predicted disk would evade detection on their SBPs if its cen- 
tral surface brightness /ve,o is fainter than 27.1 mag/rj" and its 
exponential scale length a larger than 1075. These limits trans- 
late to an apparent magnitude of 20 mag and a Ji ^28.6 mag/rj", 
slightly fainter, but probably still consistent, within the uncer- 
tainties, with the Ji predicted by the Legrand model. Thus, previ- 
ous surface photometry could neither strictly rule out nor estab- 



lish the presence of the putative ultra-LSB disk in I Zw 18. Since 
deep HST ACS narrow band imagery is unavailable for I Zw 18, 
we can not improve on the nebular line subtraction carried out 
by P02 and push previous V and R' surface photometry to fainter 
levels. 

However, CMD studies of I Zw 18 argue against the disk 
hypothesis, as none among them has revealed a uniform and ex- 
tended population of RGB stars being cospatial with the nebu- 
lar envelope, at sharp contrast to any evolved BCD studied as 
yet. Admittedly, as none among the published CMD studies for 
I Zw 18 fully appreciates the importance of ne contamination 
(except for that by IT04) and attempts a proper assessment of 
the systematic errors this can introduce, the robustness of the 
non-detection of RGB candidates in the LSB envelope can not 
be currently evaluated, neither is it clear whether CMD stud- 
ies can ever place firm constraints in this respect. The currently 
most convincing evidence for the absence of a substantial stellar 
background beneath the nebular envelope of I Zw 1 8 is based on 
the emission-line free HST images by P02 (see their discussion). 
Note that, even if an ultra-LSB stellar disk complying with the 
photometric limits placed by P02 was present, it would likely 
not significantly alter the youth interpretation for I Zw 18: The 
absolute V magnitude of the predicted disk (-1 1.4 mag), implies 
for continuous SFR over 14 Gyr a stellar mass of ~ 5.6x 10 6 Mq. 
The V magnitude of the host galaxy of I Zw 18 (i.e. excluding 
regions NW and SE) of -13.8 mag (P02) translates by SFH2 for 
an age between 0.5 and 1 Gyr to a stellar mass of - llxlO 6 M , 
i.e. about twice that of the hypothetical ultra-LSB disk. Even in 
the extreme case of a 14 Gyr old instantaneous burst, the under- 
lying disk would only be twice as heavy as the young population 
which have formed on a much shorter time scale. Consequently, 
the conclusion that I Zw 18 has formed most of its stellar mass at 
a late cosmic epoch still holds. It would be interesting to check 
whether the combined chemical output from these two stellar 
populations is consistent with the observed gas-phase metallic- 
ity in I Zw 18. 

With regard to the proposed generalization of the 'slowly 
cooking' scenario, it should be noted that the Legrand model 
faces difficulties in reproducing basic photometric properties of 
BCDs. Specifically, the host galaxy of typical BCDs is incom- 
patible with an ultra-LSB disk and, to the contrary, its central 
stellar density (>1 M pc~ 3 ) exceeds by an order of magnitude 
that of dwarf irregulars or dwarf spheroidals (P96b). Similarly, 
the mean surface brightness of the BCD host is <24 B mag/rj" 
(e.g. P96b, Amorin et al., 2009), roughly 4 mag brighter than 
the model prediction. Additionally, as pointed out in Sect. 4.1.1, 
continuous star formation models are inconsistent with the red 
host galaxy colors of most nearby BCDs. In summary, whereas 
appealing from the chemical point of view, the 'slowly cook- 
ing' scenario does not seem to be reconcilable with the struc- 
tural properties of BCDs. This scenario might, however, pro- 
vide a good approximation to quiescent late-type dwarfs, such as 
e.g. blue low surface-brightness galaxies (Ronnback & Bergvall, 
1994). 

4.3. I Zw 18 as local morphological template for rapidly 
assembling galaxies at high redshift 

The integral photometric properties of the overwhelming ma- 
jority of SF galaxies in the local universe are barely affected 
by ne. Typically, the EW(Hff+[N II]) ranges between a few ten 
A (e.g. Moustakas & Kennicutt, 2006; Koopmann & Kenney, 
2006) in normal late-type galaxies to $ 10 2 A for the majority of 
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local SFDGs (see e.g. Lee et al., 2007; Sanchez Almeida et al., 
2008). Even for BCDs, ne has a noticeable photometric impact 
(EW(Ha+[Nn])^ 200 - 500 A; Terlevich et al., 1991; Cairos 
et al., 2002; Bergvall & Ostlin, 2002; Gil de Paz et al., 2003; 
Guseva et al., 2009, among others) in their centrally confined SF 
component only (i.e. for R*^R$ F , with R$ F being of the order of 
~ 1 kpc) and is practically negligible for larger radii where the 
evolved stellar LSB host entirely dominates the line-of-sight in- 
tensity (P02, Knollmann, 2004). More generally, notwithstand- 
ing the fact that ne in local SFDGs may protrude beyond R$ F 
(Hunter & Gallagher, 1985; Gallagher et al., 1989; Hunter & 
Gallagher, 1990, 1992; Meurer et al., 1992; Marlowe et al., 
1995; Ferguson et al., 1996; Papaderos & Fricke, 1998; Martin, 
1998; Bomans, 2002; Cairos et al., 2002, among others), it has 
due to its extremely low surface brightness no impact on surface 
photometry (P02, Knollmann, 2004). As an example, even for 
BCDs with strong ongoing SF activity, corrections of integral 
photometric quantities for ne do not exceed «0.1 mag (Salzer et 
al., 1989). 

In the nearby universe, the only cases of SFDGs with ex- 
treme ne contamination are documented in a few XBCDs. 
Intense and almost galaxy-wide SF activity in these rare young 
galaxy candidates, in combination with the low surface density 
of their underlying stellar host, boost [Olll]/I5007 and Ha emis- 
sion line EWs to values of up to 2 x 10 3 A and <; 1.6 x 10 3 A, 
respectively (e.g. P98, Izotov et al, 1997a; Guseva et al., 2004; 
Izotov et al., 2006; Papaderos et al., 2006b, 2008, among oth- 
ers), i.e. of the order of the effective width of broad band fil- 
ters. Similar, though less extreme, cases are some of the recently 
discovered ultra-compact starbursting dwarfs in galaxy clusters 
(EW(Ha)~ 10 3 A, Reverte et al., 2007). 

As shown in P02, the radial Ha intensity profiles of 
BCDs comprise two characteristic components: a higher-surface 
brightness core with R*-R SF and an outer, roughly exponential 
LSB envelope with au a in the range 0.1-1 kpc. Judging from its 
aua (210 - 270 pc) and radial extent (R*~2.6 kpc), the ne halo 
of I Zw 18 is by no means exceptional among BCDs/SFDGs. 
However, I Zw 18 strikingly differs from any SFDG studied as 
yet by the fact that the galaxy itself (i.e. its stellar host, cf. Fig. 
8) is several times more compact than the nebular halo: ne dom- 
inates already for R*^6" (=3 R e ff) and reaches as far out as ~16 
stellar exponential scale lengths (Sect. 3.2). This, so far unique, 
case in the nearby universe is schematically illustrated in Fig. 
12 where the radial distribution of stars and ne in I Zw 18 is 
compared with that of normal BCDs. 

For the forthcoming discussion it is of special importance 
to recall that the surface brightness level at which in I Zw 18 
ne dominates is quite high (fi ss 23.5 mag/D"; cf Fig. 5), i.e. 
comparable to the central surface brightness of dwarf irregulars 
(e.g. Patterson & Thuan, 1996; van Zee, 2000) and dwarf ellipti- 
cals (Binggeli & Cameron, 1991, 1993), and by at least one mag 
brighter than that of Local Group dwarf spheroidals (cf. Mateo, 
1998). The ne envelope of I Zw 18 is therefore not to be con- 
fused with the extraordinarily diffuse ne in typical SFDGs. Of 
special relevance is as well the fact that the exponential ne en- 
velope contributes at least 1/3 of the total R band luminosity of 
IZw 18(P02and Sect. 3.2). 

As we shall argue next, the morphological properties of 
I Zw 18, whereas unique among nearby SFDGs, are likely typ- 
ical for rapidly assembling galaxies in the distant universe. Just 
like I Zw 18, these systems are building up their stellar compo- 
nent through starbursts or prolonged phases of strongly elevated 
specific SFR (SSFR), translating into short (a few 100 Myr) stel- 
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Fig. 12. Schematic comparison of typical Blue Compact Dwarf (BCD) 
galaxies (left) with I Zw 18 (right) with respect to the morphology and 
radial intensity distribution (upper and lower figure, respectively) of 
their stellar and nebular emission. In typical BCDs (and star-forming 
dwarf galaxies in general) the contribution of nebular emission to the 
line-of-sight intensity is practically negligible, especially outside their 
centrally confined SF component (i.e. for R*>Rsf). By contrast, in 
I Zw 18 nebular emission dominates already in the inner part of the 
galaxy (R*^3R c ff, corresponding to a surface brightness level of 23.5 
mag/D") and produces its exponential lower-surface brightness enve- 
lope that reaches as far out as ~16 stellar exponential scale lengths. 



lar mass doubling times. The cumulative output of energy and 
momentum from stellar winds and SNe during such dominant 
phases of galaxy evolution will inevitably lead to a large-scale 
gas thermalization and acceleration, with super-shells protruding 
much beyond the galaxy itself. Extended nebular halos encom- 
passing the still compact stellar component of high-SSFR proto- 
galactic systems may thus be ubiquitous in the early universe. 
The photometric structure of these galaxies could then closely 
resemble the right panel of Fig. 12, comprising a HSB core to 
within which the stellar component is confined and dominates 
and a much larger, nearly exponential, nebular LSB envelope. 
Despite cosmological dimming, the latter should be readily ac- 
cessible to observations out to z ^ 2, given that deep HST imag- 
ing and image stacking now permit studies of galaxies down to 
rest-frame surface brightnesses of ^ 26.5 - 28.5 mag/D" at those 
redshifts (e.g. Stockton et al., 2008; van Dokkum et al., 2010; 
Noeske et al., 2006). Morphological analogs to I Zw 18 may as 
well exist among compact high-SSFR galaxies at intermediate 
redshift (0.1$ z $0.8), such as e.g. Compact Narrow Emission- 
Line Galaxies (CNLEGs, Koo et al., 1994; Guzman et al., 1998), 
Luminous Compact Blue Galaxies (Guzman et al., 2003; Puech 
et al., 2006) and Green Pea (GP) galaxies (Cardamone et al., 
2009; Amorin et al., 2010; Izotov et al., 2011). Note, that the 
rest-frame EWs of GPs can in some cases almost compete with 
those of nearby XBCDs (see e.g., Izotov et al., 2001b; Papaderos 
et al., 2006b). 

Such considerations motivate a closer examination of obser- 
vational biases that the spatial segregation between stellar and 
nebular emission may introduce into studies of distant, poorly 
resolved morphological analogs of I Zw 18. By appreciation of 
the differing spatial distribution of stellar and nebular emission 
the discussion here goes beyond the framework of state-of-art 
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modeling studies (e.g. Huchra, 1977; Bergvall, 1985; Salzer et 
al., 1989; Olofsson, 1989; Kriiger et al., 1995; Fioc & Rocca- 
Volmerange, 1997; Izotov et al, 1997a; Guseva et al., 2001; Moy 
et al., 2001; Anders et al., 2004; Zackrisson et al., 2001; Panuzzo 
et al., 2003; Zackrisson et al., 2008; Kotulla et al., 2009; Martm- 
Manjon et al., 2010; Shaerer & de Barros, 2009; Finlator et al., 
201 1; Ono et al., 2010, among others) exploring in detail the im- 
pact of ne on the integral SED of SF galaxies. The predictions 
from such zero-dimensional models are to be treated with some 
caution when compared with spatially resolved observables (e.g. 
radial EW and color profiles) in order to place constraints on 
the formation history of galaxies. The minimum prerequisite for 
this approach to be valid is that nebular emission is cospatial 
with the local ionizing and non-ionizing stellar background, or, 
equivalently, that the ionizing Lyman continuum budget is re- 
processed into nebular emission on the spot. That this idealized 
picture is not invariably valid and, actually, a strong spatial anti- 
correlation between emission line EWs and stellar surface den- 
sity may evolve over time both on small and large scales has been 
shown for several XBCDs (e.g. P98; P02; Papaderos et al., 1999; 
Izotov et al., 2001a; Guseva et al., 2001; Fricke et al., 2001). 
In these cases, subtraction of synthetic ne SEDs from observed 
spectra, as done in P98, offers the only viable approach for iso- 
lating and age-dating the underlying stellar SED. 

A discussion of aperture effects on the luminosity-weighted 
SEDs of distant morphological analogs to I Zw 18 is beyond the 
scope of this study. Here, we will only focus on potential photo- 
metric biases. A first one, already described in P02, arises from 
the fact that the nebular envelope may, due to its exponentiality 
and reddish (B-R and V-R) rest frame colors, be taken as evi- 
dence for an evolved stellar disk. Automated surface photometry 
analyses of large extragalactic probes, if not including a robust 
discriminating criterion between stellar and nebular emission, 
could therefore be biased towards an enhanced frequency of 
galactic disks already at place at an early cosmic epoch. The ar- 
guably most convincing argument for this concern derives from 
I Zw 18 itself: if for one of the best studied galaxies in the local 
universe, it took three decades to realize that its exponential LSB 
envelope is not due to a stellar disk but due to ne, one may right- 
fully doubt that this would be immediately apparent in studies of 
its poorly resolved morphological analogs at high z. 

Another potential bias is owing to the fact that the superpo- 
sition of two exponential profiles of differing po and a - one 
representing the steeper star-dominated core and the other the 
shallower nebular envelope - can closely approximate a gen- 
uine Sersic profile with a high shape parameter (2< r\ ^ 5), 
thereby mimicking the SBP of a massive spheroid. The best- 
fitting Sersic exponent r\ for such a composite SBP depends both 
on the properties of its constituent exponential profiles and the 
limiting surface brightness pn m down to which Sersic models 
are fitted. As an example, the best-fitting r\ for I Zw 18 increases 
monotonically from 1.1 for pn m = 23.5 mag/D" (at the surface 
brightness where ne takes over) to 2.2 for pn m = 28 mag/D". 
This is illustrated in Fig. 13, where the observed SBP of I Zw 18 
(filled circles) is overlaid with the superposition of two exponen- 
tial components, approximating the core (open circles; profile la- 
beled core) and the ne envelope (dotted line crossing the abscissa 
at ~22 mag/D"; envl). Their sum (core+envl; orange solid-line 
curve) is then fitted by Sersic models down to a progressively 
faint p\\ m . It can be seen from Fig. 14 that the best-fitting r\ 
(squares) for the synthetic SBP is doubled when fi\ lm decreases 
from 23.5 to >27.5. The same trend is recovered when Sersic 
models are fitted directly to the observed SBP (thick curve). 
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Fig. 13. upper panel: Comparison of the R SBP of I Zw 18 (filled cir- 
cles), with a synthetic SBP (labeled core+envl) that is due to the su- 
perposition of two exponential components of differing central surface 
brightness fi and exponential scale length a. The first one, labeled core 
(open circles) approximates the stellar component that is confined to 
and dominates within the inner (R*<6") HSB part of the observed SBP. 
The second component (envl; dotted line intersecting the abscisca at 
[l ~ 22 mag/n") is a linear fit to the nebular LSB envelope. A Sersic 
model to the composite core+envl SBP (orange solid-line curve) yields 
a Sersic exponent r] «2, close to the best-fitting value of rj «2.2 for 
the observed SBP. The synthetic SBP labeled core+env2 (red solid-line 
curve) is computed by superposing on the core an exponential nebular 
envelope of equal luminosity but twice as large a as the observed enve- 
lope envl. A Sersic fit to core+env2 (thin curve) yields an rj a 5. lower 
panel: Residuals between the synthetic SBPs core+envl (thin curve) 
and core+env2 (thick curve) and their Sersic fits when the latter are 
computed down to a surface brightness level yUij m of 26.5 and 29 mag/D" 
(solid and dotted curve, respectively). It can be seen that residuals do not 
exceed 0.2 mag and 0.5 mag, respectively, i.e. they are of the order of 
or smaller than typical uncertainties in SBPs of intermediate-to-high z 
galaxies at rest frame surface brightness levels /j S26.5 mag/D". 



It is worth checking how the Sersic 77 vs. p e g relation may 
change in the case of an equally luminous but shallower nebular 
envelope. For this, we superpose to the core a component (env2) 
whose po and a is by, respectively, 1 .5 mag fainter and a factor 
of 2 larger than that in the observed component envl. The Sersic 
fit to the core+env2 profile (red solid-line curve) is included in 
Fig. 13 with the thin-dark curve. As apparent from Fig. 14, in 
this case rj shows a steeper dependence on pn m , increasing to 
~ 4 already for p\\ m 26.5 mag/D" and leveling off to ~5 at 
fainter levels. 

Consequently, Sersic fits can readily lead to the misclassi- 
fication of a compact high-SSFR galaxy as a massive elliptical 
and, quite counter-intuitively, deeper photometry exacerbates the 
problem. Note that the residuals between SBP and model (lower 
panel of Fig. 13) are S=0.5 mag for core+env2 and just $0.2 mag 
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Fig. 14. Best-fitting Sersic exponent r] vs /Ji; m for the SBPs in Fig. 
13. The thick-gray curve and the filled squares show, respectively, the 
variation of rj for the observed R band SBP of I Zw 1 8 and the two- 
component approximation to it, labeled core+envl. The variation of rj 
with fi e g for the synthetic SBP core+env2 is shown with filled circles. 



for core+envl, i.e. of the order of lcr uncertainties in SBPs of 
intermediate-to-high z galaxies (cf e.g. Noeske et al., 2006), i.e. 
small enough to go undetected. In practice, even when PSF con- 
volution effects are fully accounted for, the pseudo-Sersic profile 
of a compact diskless high-SSFR galaxy is barely distinguish- 
able from the Sersic profile of a massive galaxy spheroid. 

Evidently, since extended ne can drastically affect an SBP 
as a whole, it also impacts virtually all secondary photometric 
parameters that are derivable from it (e.g. the effective radius 
and Gini coefficient, and the various light concentration indices 
commonly used in galaxy quantitative morphology studies). 

Thirdly, the ne luminosity fraction in I Zw 18 (>l/3), if typ- 
ical for its higher-z analogs, translates into a systematic error of 
>0.4 mag in galaxy scaling relations involving total magnitudes. 
These range from the Tully-Fisher relation to all relations com- 
paring luminosity with e.g. metallicity, diameter, mean surface 
brightness and velocity dispersion. Moreover, errors in galaxy 
luminosity propagate, potentially amplified, in stellar mass de- 
terminations using theoretical mass-to-light ratios or SED fit- 
ting. An investigation of this issue was recently presented by 
Izotov et al. (2011): these authors has shown that masses com- 
puted from spectral fits to the SED continuum of high-SSFR 
galaxies at intermediate z can be overestimated by a factor of up 
to ~4. This is not primarily due to the luminosity contribution but 
due to the red spectral slope of the nebular continuum (see e.g. 
Kriiger et al., 1995) which drives SED fitting towards solutions 
invoking a much too high luminosity fraction from old stars. In 
galaxy assembly studies covering a wide range in z, downsizing 
effects may further complicate the aforementioned biases, since 
ne is expected to affect galaxy populations of different mass over 
different timescales. 

We next turn to the color contrast 6 ce between the HSB core 
and the nebular LSB envelope. This quantity can readily be de- 
termined from SBPs, or within concentric apertures, and pro- 
vides a handy proxy to radial color gradients in galaxies, thus 
a first classification guess. As we will show below, 6 ce offers 
for certain z intervals a powerful diagnostic for identifying high- 
SSFR galaxies with morphological properties analogous to those 
I Zw 18. Within other z intervals, however, and depending on the 
colors considered, a superficial interpretation of <S ce can further 
aggravate the above discussed galaxy misclassification biases. 



In computing 6 ce and its variation with z, we approximated 
the spectrum of I Zw 18 with synthetic stellar + ne SEDs from 
Pegase 2, referring to SFH2 and a metallicity Z=0.0004. In these 
models, the properties of the HSB core (R*<6") are well repro- 
duced by a synthetic SED for an age f=100 Myr. This yields col- 
ors of V-R=0.2, V-/ = -0.15, R-I = -0.37 mag, and an EW(Hcr) 
of 670 A, in good agreement with the observed values (Izotov 
et al., 2001a; Vflchez & Iglesias-Paramo, 1998, P02). As for 
the envelope (6"<R*^20"), we adopted the SED from the same 
model for f=0, i.e. a purely ne spectrum. Its colors (fi-V=0.28, 
V-R=QA1, B-R=0J, V-I = -0.65 and R-I = -1.1 mag) provide 
as well a good match to the data. 

The variation of the B-J, V-K, V-R, V-I and R-I colors 
of the core and the envelope as a function of z is shown in the 
upper two panels of Fig. 15. It can be seen that the envelope 
(middle panel) shows particularly large color variations, as dif- 
ferent strong emission lines shift into the transmittance window 
of various filters, depending on z. As already evident from Sect. 
3.2, a local analog to I Zw 18 is easily identifiable by its blue 
nebular envelope and large (0.8 mag) S ce (lower panel) both with 
respect to V-I and R-I. This is also the case for the redshift 
range 0.15 55 z ^ 0.3 where the envelope appears much red- 
der than the core in V-I but bluer in B-V. Other distinct peaks 
in ||5 ce || with respect to various optical or optical-NIR colors 
are apparent for e.g. z ~ 0.42, 0.55 and 0.9. A comparison of 
the upper and lower panel of Fig. 15 shows that the 6 ce exhibits 
much larger variations (||£ ce ||-1.6 mag) than the HSB core, i.e. 
it is a far more sensitive indicator of extended ne than integral 
(luminosity-weighted) colors that are primarily driven by the 
core. For example, at z ~ 0.27, the V-I and B-V colors of the 
HSB core are equal to within <0.3 mag, whereas their <S ce is dif- 
fering by up to x 1.3 mag. 

It is worth pointing out that, for certain redshift windows 
and color indices, the (5 ce may point towards diametrically dif- 
ferent views on the nature of an I Zw 18 analog. For instance, at 
0. 15 ^ z ^ 0.3 the large V-I core-to-envelope color contrast (6 ce 
~ 0.8 mag) and the moderately blue colors of the core (0.5 mag) 
superficially suggest an old disk hosting nuclear SF activity. The 
opposite conclusion could be drawn from the <5 ce in B-V («0.5 
mag) which may be taken as evidence for a very young stellar 
disk encompassing a slightly older core, in line with the inside- 
out galaxy growth interpretation. In either case, the exponential 
envelope, usually interpreted as stellar disk, would nicely aug- 
ment the erroneous conclusion. Several similar cases for various 
color combinations and redshifts can be read off Fig. 15. Of im- 
portance is also the fact that for other z's (e.g. 0.1& z 550.15, 0.3, 
0.75, 1.05) the S ce vanishes to SS0.2 mag for some colors, having 
little discriminating power between stellar and nebular emission, 
and superficially suggesting a nearly uniform stellar age. 

With regard to all concerns above, one might counter-argue 
that photometric k corrections would rectify SBPs and color pro- 
files, and eliminate pitfalls in the interpretation of 5 ce . However, 
state-of-art k corrections, mostly tailored to stellar SEDs and ap- 
plied to a morphological analog of I Zw 18 as a whole, regard- 
less of its physically distinct radial zones, would most probably 
aggravate the problem in a barely predictable manner. 

In summary, the case of I Zw 18 stands as a warning bench- 
mark to studies of high-SSFR galaxies near and far. It reminds 
us of the fact that nebular emission in SF galaxies is not always 
cospatial with the underlying local ionizing and non-ionizing 
stellar background, neither has to scale with its surface den- 
sity. Quite the contrary, in galaxies with high SSFR (i.e. forming 
young galaxies and/or systems rapidly assembling their stellar 
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Fig. 15. Variation of the observed color as a function of redshift z for a 
template galaxy with the properties of I Zw 18. The model consists of 
a high-surface brightness core to within which stellar emission is con- 
fined and dominates (upper panel), and a larger lower-surface bright- 
ness nebular envelope (middle panel). The colors of the core and the 
envelope yield for z = a good match to the observed values for the 
core (R*<6") and the nebular envelope (6"<1?*S20") of I Zw 18 (cf 
Sect. 3.2). The lower panel shows the dependence of the color contrast 
<5 CC between the core and the envelope as a function of redshift. 



mass during dominant phases of their evolution) nebular emis- 
sion is plausibly expected to extend far beyond the galaxy it- 
self and to significantly contribute to its total luminosity. This, 
in connection with the empirical evidence that nebular emission 
forms an exponential envelope, conspires in potentially impor- 
tant observational biases in studies of moderately resolved mor- 
phological analogs to I Zw 18 which are likely ubiquitous in the 
early universe. 

Interesting in this context is also the detection of exponen- 
tial Lya halos around low (Hayes et al., 2007) and high-z star- 
forming galaxies (Steidel et al., 2011), resulting in SBPs strik- 
ingly similar to those of I Zw 18. 

In principle, any strong thermal or non-thermal central 
source of energy and momentum could generate a large nebular 
envelope and a pseudo bulge-disk luminosity profile. The fore- 
going discussion is therefore of relevance to a wider range of 
topics, as e.g. the co-evolution of AGNs with their host galaxies 
or the properties of ionized halos around powerful radio galaxies 
(e.g., Villar-Martin et al., 2003, and references therein). 

The prospect of spatially resolved studies of galaxy forma- 
tion in the faraway universe with next-generation observing fa- 



cilities calls for theoretical guidance on the properties and time 
evolution of nebular halos around protogalactic systems. Some 
questions (see also P02) to computational astrophysics include: 
i) in which way are the photometric properties of the nebular 
envelope (e.g. its po and a), and their temporal evolution, re- 
lated to the specifics of energy production (e.g. the SFH and 
the ionizing SED), the initial geometry and kinematics of the 
protogalactic gas reservoir and the physical conditions in the 
intergalactic medium (e.g. external pressure, ambient ionizing 
field from multiple mutually 'illuminating' protogalactic units)? 
More specifically, ii) does (for a given set of environmental con- 
ditions) the core-envelope radial intensity pattern scale in a ho- 
mologous manner, i.e. is the a of the nebular envelope invari- 
ant? If not, iii) what does the core-to-envelope luminosity and 
normalized Rs^/a ratio tell us about e.g. the recent SFH and cur- 
rent SSFR? Additionally, iv) how is the photometric, chemical 
and kinematical evolution of the nebular envelope related to the 
escape probability of Ly continuum and Lya photons? 

5. Summary and conclusions 

We used the entire set of archival HST ACS broad band imaging 
data for I Zw 1 8 and its fainter component I Zw 1 8 C to study the 
photometric structure of this nearby (19 Mpc) dwarf galaxy pair 
to unprecedently faint surface brightness levels (p ^29 mag/D"). 

The main results from this study may be summarized as follows: 

i) Radial color profiles reveal very blue and practically constant 
colors (0±0.05 mag) for I Zw 18 C down to p ~ 27.6 mag/D", 
and a previously undisclosed, slightly redder (V-/=0.2+0.08 
mag, R-I=Q. 16+0.09 mag) stellar component in its extreme pe- 
riphery (27.6 - 29 mag/D"). We have verified that these blue 
colors do not merely reflect the luminosity-weighted average of 
blue young stellar clusters with a faint red stellar background 
but that they are characteristic for the unresolved host galaxy 
of I Zw 18 C. We argue that the buildup of the photometri- 
cally dominant stellar component of I Zw 18 C has occurred in 
a largely sequential mode, through a star-forming (SF) process 
that has likely started r ~100 Myr ago at the redder northwest- 
ern tip of the galaxy and propagated with a mean velocity of ~20 
km s _1 to its bluer southeastern tip. 

The photometric properties of the extreme periphery of 
I Zw 18 C are entirely consistent with this formation scenario, 
if the effect of stellar diffusion is taken into account. Radial mi- 
gration of newly forming stars with a mean velocity of «4 km 
s over ~ r, and the associated stellar mass filtering effect de- 
scribed in Papaderos et al. (2002, hereafter P02), can naturally 
account for the slightly redder colors, absence of nebular emis- 
sion (ne) and topological properties of the stellar outskirts of 
I Zw 18C. A faint ancient stellar substrate can not be ruled out, 
even though our analysis does not lend observational support to 
its presence. 

ii) In I Zw 18 (i.e. the main body) severe contamination of broad 
band colors by ne prevents a conclusive age dating of the stel- 
lar component almost everywhere. Therefore, even though our 
combined images are the deepest presently available, we can 
not improve on the age-dating analysis by P02 who have two- 
dimensionally subtracted strong nebular emission lines from 
broad band HST WFPC2 images to isolate and age-date the 
residual underlying stellar background. However, the colors of 
the reddest quartile of region a>, a region at the southeastern tip 
of I Zw 1 8 with comparatively weak ne, were determined to be 
in good agreement with those previously inferred by P02. These 
colors, if representative for the host galaxy of I Zw 18, imply on 
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the basis of a continuous or exponentially decreasing star forma- 
tion history (SFH) that I Zw 18 has formed most of its stellar 
mass at a late cosmic epoch. This, together with our conclusions 
under i), supports the view that both I Zw 18 and I Zw 18 C are 
cosmologically young objects that have undergone a nearly syn- 
chronous evolution. 

iii) We show that ~2Q% of the isophotal area of I Zw 18 at 25 
mag/D" is severely affected by ne, thus inaccessible to age dat- 
ing studies via broad band colors and color-magnitude diagrams 
(CMDs). The local impact of ne manifests itself i.a. in a combi- 
nation of red (~0.5 mag) V-R with extremely blue (-0.4 . . . -0.8 
mag) V-I and R-I colors. Nebular emission shows considerable 
sub-structure, with numerous clumps and overlapping shells, and 
little spatial correlation with the local stellar background. It thus 
can not be treated as a uniform foreground emitting layer and ac- 
curately subtracted out using standard point source photometry 
algorithms. This likely results in substantial random and system- 
atic errors that might not be fully accounted for by the standard 
CMD error budget. Further potential sources of systematic un- 
certainties stem from spatial displacements (50 - 100 pc) be- 
tween the intensity maxima of the [Olll]/i/l4959,5007 and Ha 
emission lines along ionized gas shells. These may differentially 
affect the local background determination in various filters, caus- 
ing an artificial reddening of CMD point sources in the interior 
of nebular shells (and vice versa). 

iv) Based on the extraordinarily deep combined HST ACS im- 
ages, we have been able to study the exponential low-surface 
brightness (LSB) envelope of I Zw 18 out to its extreme pe- 
riphery using both surface brightness profiles (SBPs) and color 
maps. These reveal uniform colors of V-/?«0.55 mag, V-I~-l 
mag and R-Ix-l .4 mag all over the LSB component of I Zw 18, 
corroborating the previous conclusion by P02 that this luminos- 
ity envelope is not due to a stellar disk but due to extended ne. 
Specifically, our analysis indicates that ne dominates the line- 
of-sight intensity beyond 3 effective radii (or, equivalently, for 
fi >23.5 mag/D") and extends as far out as ~16 stellar exponen- 
tial scale lengths. The overall picture emerging from our anal- 
ysis is therefore that I Zw 18 is a cosmologically young object 
that consists of a compact, high-surface brightness (HSB) core, 
within which stellar emission is confined and dominates, and a 
much larger nebular LSB envelope. 

v) We argue that the morphological properties of I Zw 18, while 
unique among nearby SF dwarf galaxies, are probably typical 
among distant young galaxies in dominant phases of their evo- 
lution, during which they assemble their stellar component at 
high specific star formation rates (SSFRs). The prodigious ener- 
getic output during such phases of rapid stellar mass growth is 
expected to result into a large-scale gas ionization and acceler- 
ation, with ne protruding much beyond the still compact stellar 
galaxy host. These systems could thus bear strong morpholog- 
ical resemblance to I Zw 18, comprising a compact core that 
is dominated by stellar emission and a much larger exponential 
nebular envelope. A question of considerable interest is there- 
fore, how the nebular envelope could impact photometric stud- 
ies of moderately resolved morphological analogs of I Zw 18 at 
higher z. 

A potential bias, already discussed in P02, arises from the 
fact that the nebular envelope mimicks due to its exponentiality 
and red B-R color an evolved stellar disk. Here, by using I Zw 1 8 
as a template, we extend previous considerations: 

v.l We point out that the superposition of two exponential 
components of differing central surface brightness and scale 
length, approximating the core and the envelope of a distant 



I Zw 18 analog, may be barely distinguishable from a genuine 
Sersic profile with an exponent 2% rj 5=5. Therefore, Sersic mod- 
els offer, in the specific context, a poor diagnostic for disen- 
tangling compact high-SSFR protogalaxies from massive galaxy 
spheroids. 

v.2 Nebular emission contributes at least 1/3 of the total lu- 
minosity of I Zw 18 (P02 and this study). This luminosity frac- 
tion, if typical for its higher-z analogs, translates into a system- 
atic error of ^0.4 mag in all galaxy scaling relations involving 
luminosities (e.g., the Tully-Fisher relation, and relations be- 
tween luminosity and metallicity, diameter, velocity dispersion). 
Evidently, errors in total luminosities propagate into errors in 
galaxy mass determinations via theoretical mass-to-light ratios 
or SED fitting techniques. Moreover, since extended ne can dras- 
tically affect a galaxy SBP as a whole, it also affects practically 
all secondary photometric quantities that are derivable from it 
(e.g. the effective radius or various light concentration indices 
used for quantitative galaxy morphology studies). 

v.3 We investigate the variation of the color contrast <5 ce be- 
tween the star-dominated core and the surrounding nebular enve- 
lope as a function of z. This task is motivated by the fact that 6 ce 
provides a handy proxy to radial color gradients in galaxies and 
a valuable galaxy classification tool. We show that for certain z 
intervals, this quantity offers a powerful diagnostic for the iden- 
tification of moderately resolved I Zw 18 analogs. Within other 
z intervals, however, and depending on the color indices con- 
sidered, a superficial interpretation of <5 ce can further enhance 
galaxy misclassification biases stemming from SBP fitting (cf 
v.l) and potentially impact our understanding of galaxy assem- 
bly over time. State-of-art k corrections applied to distant mor- 
phological analogs to I Zw 18 as a whole, i.e. regardless of their 
physically distinct radial zones, may aggravate observational and 
interpretation biases in a barely predictable manner. 

In the era of spatially resolved studies of galaxy formation 
in the early universe with next-generation observing facilities, 
a better theoretical understanding of the rise and fall of nebular 
galaxy halos over cosmic time appears to be crucially important. 
In this respect, some questions to computational astrophysics are 
formulated. 
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